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Abstract 

This  paper  presents  the  test  results  and  life  modeling  of  special  calendar-  and  cycle-life  tests  conducted  on  18650-size  generation  1  (Gen  1) 
lithium-ion  battery  cells  (nominal  capacity  of  0.9  Ah;  3.0-4. 1  V  rating)  developed  to  establish  a  baseline  chemistry  and  performance  for  the 
Department  of  Energy  sponsored  advanced  technology  development  (ATD)  program.  Electrical  performance  testing  was  conducted  at  the 
Argonne  National  Laboratory  (ANL),  Sandia  National  Laboratory  (SNL)  and  the  Idaho  National  Engineering  and  Environmental  Laboratory 
(INEEL). 

As  part  of  the  electrical  performance  testing,  a  new  calendar-life  test  protocol  was  used.  The  test  consisted  of  a  once  per  day  discharge  and 
charge  pulse  designed  to  have  minimal  impact  on  the  cell  yet  establish  its  performance  over  a  period  of  time  such  that  the  calendar-life  of  the 
cell  could  be  determined.  The  calendar-life  test  matrix  included  two  states-of-charge  (SOCs)  (i.e.  60  and  80%)  and  four  test  temperatures  (40, 
50,  60  and  70  °C).  Discharge  and  regen  resistances  were  calculated  from  the  test  data.  Results  indicate  that  both  the  discharge  and  regen 
resistances  increased  non-linearly  as  a  function  of  the  test  time.  The  magnitude  of  the  resistances  depended  on  the  temperature  and  SOC  at 
which  the  test  was  conducted.  Both  resistances  had  a  non-linear  increase  with  respect  to  time  at  test  temperature.  The  discharge  resistances  are 
greater  than  the  regen  resistances  at  all  of  the  test  temperatures  of  40,  50,  60  and  70  °C.  For  both  the  discharge  and  regen  resistances,  generally 
the  higher  the  test  temperature,  the  lower  the  resistance. 

The  measured  resistances  were  then  used  to  develop  an  empirical  model  that  was  used  to  predict  the  calendar-life  of  the  cells.  This  model 
accounted  for  the  time,  temperature  and  SOC  of  the  batteries  during  the  calendar-life  test.  The  functional  form  of  the  model  is  given  by: 
R(t,  T,  SOC)  =  A(T ]  SOC)/7/)  +  B(T ]  SOC),  where  t  is  the  time  at  test  temperature,  T  the  test  temperature  and  SOC  the  SOC  of  the  cell  at  the 
start  of  the  test.  A(T,  SOC)  and  B(T,  SOC)  are  assumed  to  be  functions  of  the  temperature  and  SOC;  F  is  assumed  to  only  be  a  function  of  the 
time  at  test  temperature.  Using  curve-fitting  techniques  for  a  number  of  time-dependent  functions,  it  was  found  that  both  the  discharge  and 
regen  resistances  were  best  correlated  with  F(t )  having  a  square-root  of  test  time  dependence.  These  results  led  to  the  relationship  for  the 
discharge  and  regen  resistances:  R(t ,  T,  SOC)  =  A(T,  SOC)?1/2  +  B(T,  SOC).  The  square-root  of  time  dependence  can  be  accounted  for  by 
either  a  one-dimensional  diffusion  type  of  mechanism,  presumably  of  the  lithium-ions  or  by  a  parabolic  growth  mechanism  for  the  growth 
of  a  thin-film  solid  electrolyte  interface  (SEI)  layer  on  the  anode  and/or  cathode.  The  temperature  dependence  of  the  resistance  was 
then  investigated  using  various  model  fits  to  the  functions  A(T,  SOC)  and  B(T,  SOC).  The  results  of  this  exercise  lead  to  a  functional  form 
for  the  temperature  dependence  of  the  fitting  functions  having  an  Arrhenius-like  form:  A(T,  SOC)  =  a(SOC){exp[Z?(SOC)/7]}  and 
R(T,SOC)  =  c(SOC){exp[<i(SOC)/r]},  where  a  and  c  are  constants,  and  b  and  d  are  related  to  activation  energy  (Eh  and  Ed)  by  using 
the  gas  constant  (. R )  such  that  b  =  Eb/R  and  d  =  Ed/R.  The  functional  form,  therefore,  for  the  discharge  and  regen  resistances,  including  the 
SOC,  is  then:  R(t,  T,  SOC)  =  a(SOC){exp[6(SOC)/7]}A2  +  c(SOC){exp[<7(SOC)/7]}.  The  a,  b ,  c  and  d  parameters  are  explicitly  shown  as 
being  functions  of  the  SOC.  However,  due  to  the  lack  of  testing  at  SOC  values  other  than  60  and  80%  SOC,  the  exact  form  of  the  SOC  dependence 
could  not  be  determined  from  the  experimental  data.  The  values  of  a,  b,  c  and  d  were  determined,  thus  permitting  the  function  R(t,  T. \  SOC)  to  be 
used  to  correlate  the  discharge  and  regen  data  and  to  predict  what  the  resistances  would  be  at  different  test  times  and  temperatures. 

This  paper  also  presents,  discusses  and  models  the  results  of  a  special  cycle-life  test  conducted  for  a  period  of  time  at  specified  temperatures 
of  40,  50,  60  and  70  °C.  This  test,  consisting  of  specified  discharge  and  charge  protocols,  was  designed  to  establish  the  cycle-life  performance 
of  the  cells  over  a  time  interval  such  that  their  cycle-life  could  be  determined.  The  cycle-life  test  was  conducted  at  60%  SOC,  with  SOC  swings 
of  A3,  A6  and  A9%.  During  the  cycle-life  test,  the  discharge  and  regen  resistances  were  determined  after  every  100  test  cycles.  The  results  of 
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the  cycle-life  testing  indicate  that  both  the  discharge  and  regen  resistances  increased  non-linearly  as  a  function  of  the  test  time  at  each  A% 
SOC  test.  The  magnitude  of  the  resistances  and  the  rate  at  which  they  changed  depended  on  the  temperature  and  A%  SOC  at  which  the  test  was 
conducted.  Both  resistances  had  a  non-linear  increase  with  respect  to  time  at  test  temperature,  i.e.  as  the  number  of  test  cycles  increased  the 
discharge  and  regen  resistances  increased  also.  For  a  given  A%  SOC  test,  the  discharge  resistances  are  greater  than  the  regen  resistances  at  all 
of  the  test  temperatures  of  40,  50,  60  and  70  °C.  For  both  the  discharge  and  regen  resistances,  generally  the  higher  the  test  temperature,  the 
lower  the  resistance.  At  each  of  the  four  test  temperatures,  the  magnitude  of  the  discharge  and  regen  resistances  was  generally  in  the  following 
order:  A3%  SOC  >  A9%  SOC  >  A6%  SOC,  but  the  ordering  was  dependent  on  test  time. 

A  model  was  also  developed  to  account  for  the  time,  temperature,  SOC  and  A%  SOC  of  the  batteries  during  the  cycle-life  test.  The 
functional  form  of  the  model  is  given  by  R(t ,  T,  SOC,  A%  SOC)  =  A(T,  SOC,  A%  SOC)F(t)  +  B(T,  SOC,  A%  SOC)  where  t  is  the  time  at  test 
temperature,  T  the  test  temperature,  SOC  the  SOC  of  the  cell  at  the  start  and  end  of  the  test  and  A%  SOC  the  SOC  swing  during  the  test.  A(T, 
SOC,  A%  SOC)  and  B(T,  SOC,  A%  SOC)  are  assumed  to  be  functions  of  the  test  temperature,  SOC  and  A%  SOC  swing.  F(t)  is  assumed  to 
only  be  a  function  of  the  test  time  at  test  temperature.  Using  curve-fitting  techniques  for  a  number  of  time-dependent  functions,  it  was  found 
that  both  the  discharge  and  regen  resistances  were  best  correlated  by  a  square-root  of  test  time  dependence.  These  results  led  to  the  relationship 
for  the  discharge  and  regen  resistances  having  the  form  R(t ,  T,  SOC,  A%  SOC)  =  A(T,  SOC,  A%  SOC)f'/2  +  B(T,  SOC,  A%  SOC).  This 
model  is  essentially  the  same  as  used  to  analyze  the  calendar-life  test  data.  The  temperature  dependence  of  the  resistance  was  then  investigated 
using  various  model  fits  to  the  functions  A(T)  and  B(T).  The  results  of  this  exercise  lead  to  a  functional  form  for  the  functions  having  again  an 
Arrhenius-like  form:  A(T )  =  a[exp(b/T)]  and  B(T )  =  c[exp(d/T)]  where  a  and  c  are  constants,  and  b  and  d  are  related  to  activation  energies. 
The  functional  form,  therefore,  for  the  discharge  and  regen  resistances  including  the  SOC  and  A%  SOC  is  R(t,  T ,  SOC,  A%  SOC)  =  a(SOC, 
A%  SOC){exp[7?(SOC,  A%  SOC)/7]  R1/2  +  c(SOC,  A%  SOC){expfi/(SOC,  A%  SOC)/7] }.  The  a,  b,  c  and  d parameters  are  explicitly  shown 
as  being  functions  of  the  SOC  and  the  A%  SOC.  However,  due  to  the  lack  of  testing  at  SOC  values  other  than  60%  SOC,  the  exact  form  of  the 
SOC  dependence  could  not  be  determined  from  the  experimental  data.  In  addition,  no  model  was  found  that  consistently  correlated  the 
observed  resistance  changes  with  the  A%  SOC  of  the  tests.  Eliminating  the  SOC  and  A%  SOC  from  the  resistance  function,  the  function  R(t, 
T)  was  then  used  to  correlate  the  discharge  and  regen  resistances  data.  This  model  also  allows  the  prediction  of  what  the  resistances  would  be 
at  different  test  times  at  a  particular  A%  SOC  test  condition  and  temperature. 
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1.  Introduction 

The  DOE  Office  of  Advanced  Automotive  Technologies 
(OAAT),  through  the  partnership  for  a  new  generation  of 
vehicles  (PNGV)  advanced  technology  development  (ATD) 
program,  is  engaged  in  the  study  of  18650-size  lithium-ion 
cells  to  investigate  the  physical/chemical  mechanisms  caus¬ 
ing  performance  degradation  over  the  life  of  these  batteries. 
Concentrating  on  high-power  battery  development,  the  ATD 
program  supports  the  PNGV,  a  government-industry  partner¬ 
ship  striving  to  develop,  by  2004,  a  mid-size  passenger 
vehicle  capable  of  achieving  up  to  three  times  the  fuel 
economy  of  today’s  vehicles,  while  adhering  to  future 
emissions  standards  and  maintaining  affordability,  perfor¬ 
mance,  safety  and  comfort.  The  ATD  program  addresses 
these  technical  challenges  through  five  major  program  areas: 
(1)  baseline  cell  development;  (2)  diagnostic  evaluations;  (3) 
electrochemical  improvements;  (4)  advanced  materials 
development;  and  (5)  low-cost  packaging.  The  major  objec¬ 
tive  of  the  work  is  to  determine  the  causes  of  power  fade 
after  the  ATD  program  designed  cells  are  exposed  to  ele¬ 
vated  temperatures  and  tested  under  various  electrical  per¬ 
formance  evaluation  tests.  Another  objective  is  to  develop 
diagnostic  analysis  methods  that  can  be  used  to  determine 
the  physical/chemical  causes  for  cell  degradation. 

This  paper  presents  the  electrical  performance  of  the 
generation  1  (Gen  1)  lithium-ion  cells  developed  by  the 
ATD  program  [1]  during  special  calendar-  and  cycle-life 
testing  conducted  at  various  temperatures  and  state-of- 


charge  (SOC)  [2].  All  tests  were  conducted  either  in  the 
Energy  Storage  Testing  (EST)  Laboratory,  which  is  part  of 
the  Transportation  Technologies  and  Infrastructure  Depart¬ 
ment  at  the  Idaho  National  Engineering  and  Environmental 
Laboratory  (INEEL),  at  the  Lithium  Battery  R&D  Depart¬ 
ment  1521  at  Sandia  National  Laboratory  (SNL)  and  at  the 
electrochemical  energy  technology  program  at  Argonne 
National  Laboratory  (ANL). 

One  of  the  main  points  of  this  paper  is  to  present,  discuss 
and  model  the  calendar-life  test  data  on  the  lithium-ion 
batteries  developed  by  the  ATD  program.  Calendar-life  is 
of  great  importance  as  one  of  the  PNGV  goals  is  to  develop 
battery  systems  having  a  lifetime  of  15  years  if  they  are  to  be 
a  viable  energy  storage/source  for  the  next  generation  of 
vehicles.  Calendar-life  has  been  defined  by  two  general 
statements.  The  US  ABC  Electric  Vehicles  Battery  Test 
Procedures  Manual,  Revision  2  [3]  defines  calendar-life 
as  “the  length  of  time  a  battery  can  undergo  some  defined 
operation  before  failing  to  meet  its  specified  end-of-life 
criteria.”  This  definition  is  rather  nebulous.  The  specific 
definition  of  some  defined  operation  would  need  to  be 
clarified  and  stipulated  in  the  test  plan  for  a  given  battery 
system.  The  PNGV  Battery  Test  Manual,  Revision  2  [4],  is 
also  rather  terse  in  its  definition  of  calendar-life:  “this  test  is 
designed  to  permit  the  evaluation  of  cell  degradation  as  a 
result  of  the  passage  of  time  with  minimal  usage.  It  is  not  a 
pure  shelf-life  test,  because  the  cells  under  test  are  periodi¬ 
cally  subjected  to  reference  discharges  to  determine  the 
changes  (if  any)  in  their  performance  characteristics.” 
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The  ATD  program  [1]  has  attempted  to  more  clearly  define 
what  is  meant  by  a  calendar-life  test.  The  calendar-life  test  as 
used  in  the  ATD  program  has  been  designed  to  have  a 
minimal  impact  on  the  cell,  yet  still  subject  the  cell  to  a 
well- specified  charge/discharge  test  sequence  over  a  well- 
defined  period  of  time. 

The  other  topic  of  this  paper  is  to  present,  analyze  and 
model  cycle-life  test  data  on  the  cells  developed  by  the  ATD 
program.  Cycle-life  is  also  of  great  importance,  as  battery 
systems  are  expected  to  have  a  lifetime  of  15  years.  Two 
general  statements  have  defined  cycle-life.  The  US  ABC 
Electric  Vehicles  Battery  Test  Procedures  Manual,  Revision 
2  [3]  defines  cycle-life  as  “the  number  of  cycles,  each  to  a 
specified  discharge  and  charge  termination  criteria,  such  as 
depth-of-discharge  (DOD),  under  a  specified  charge  and 
discharge  regime,  that  a  battery  can  undergo  before  failing 
to  meet  its  specified  end-of-life  criteria.”  This  definition 
clearly  depends  on  the  specifics  of  the  test  protocols  as 
specified  in  the  test  plan  for  a  given  battery  system.  The 
PNGV  Battery  Test  Manual,  Revision  2  (and  3),  discusses 
various  cycle-life  tests  that  may  be  used  when  testing 
batteries  [4] .  The  cycle-life  test  as  used  in  the  ATD  program 
has  been  designed  and  specified  as  applicable  to  the  Gen  1 
cells  developed  for  this  program. 

Discussion  of  the  terminology  used  in  this  report  can  be 
found  in  [3],  USABC  Electric  Vehicle  Battery  Test  Proce¬ 
dures  Manual,  Revision  2  and  in  [4],  PNGV  Battery  Test 
Manual,  Revision  2  (and  3).  The  entire  test  procedure  used  to 
test  the  ATD  Gen  1  cells  is  not  reproduced  herein  as  much  as 
it  is  given  in  [2],  PNGV  Test  Plan  for  ATD  18650  Gen  1 
Lithium-Ion  Cells. 

The  results  presented  in  this  paper  are  more  extensively 
presented  in  two  INEEL  reports  [5,6].  A  paper  concerning 
the  rates  of  the  area-specific  impedance  increase  and  power 
fade  resulting  from  the  accelerated  calendar-  and  cycle-life 
tests  of  the  ATD  Gen  1  cells  has  also  been  published  [7]. 

2.  Experimental 

2.7.  Description  of  ATD  program  Gen  1  lithium-ion 
batteries 

The  baseline  lithium-ion  cells  had  the  following  specifi¬ 
cations,  as  developed  by  ANL  for  the  ATD  program.  Cells 
produced  with  these  specified  materials  are  referred  too  as 
Gen  1  cells. 


Positive  electrode  (cathode)  (copper  current  collector) 
LiNio.8Coo.2O2  (Sumitomo)  (84  wt.%) 

Electronic  additive:  acetylene  black  (4  wt.%)  +  SFG-6 
graphite  (Timcal)  (4  wt.%) 

Binder:  polyvinylidene  fluoride  (-CH2-CF2-)n, 
(PVDF)  (Kureha  KF-1100)  (8  wt.%) 

Negative  electrode  (anode)  (aluminum  current  collector) 


Blend  of  MCMB-6-2800  graphite  (Osaka  gas)  (75  wt.%) 
and  SFG-6  (Timcal)  (16  wt.%) 

Binder:  PVDF  (Kureha  C)  (9  wt.%) 

Electrolyte 

1  M  LiPF6/ethylene  carbonate  (EC)  (1  wt.%)  -f  diethyl 
carbonate  (DEC)  (1  wt.%) 

Separator 

Polyethylene  (PE)  Celgard  separator  (37  pm  thick) 


Three  hundred  cells  were  assembled  (18650-size; 
64.9  mm  high,  18.12  mm  diameter)  and  shipped  to  various 
national  laboratories  [ANL,  Brookhaven  National  Labora¬ 
tory  (BNL),  INEEL,  Lawrence  Berkeley  National  Labora¬ 
tory  (LBNL)  and  SNL]  for  electrical  performance  testing 
and  physical/chemical  diagnostic  analysis.  The  cell  distri¬ 
bution  is  given  in  the  test  plan  [2] .  For  the  various  tempera¬ 
ture  tests,  controlled  temperature  chambers  having  both 
heating  and  cooling  capabilities  were  used.  Temperature 
control  was  specified  to  be  ±3  °C. 

2.2.  Electrical  performance  ratings  of  the  ATD 
program  Gen  1  lithium-ion  cells 


The  ATD  Gen  1  cell  limits  are  the  following. 


Nominal  rated  cell 

0.9  Ah  at  a  Cf  1  discharge  rate 

capacity 

Discharge 

Minimum  discharge 

3.0  V 

voltage 

Maximum  discharge 

2.0  A  continuous;  7.2  A  (8  C)  for 

current 

up  to  an  18  s  pulse;  and  13.5  A 
(15  C)  for  up  to  a  2  s  pulse 

Maximum  discharge 

temperature 

70  °C 

Charge  and  regen 

Maximum 

4.1  V  continuous;  4.3  V  for  up 

charge/regen  voltage 

to  a  2  s  pulse 

Maximum 

0.9  A  continuous  charge 

charge/regen  current 

current;  12  A  maximum  regen 
current  for  up  to  a  2  s  pulse 

Maximum  charge 

40  °C 

temperature 

Maximum  regen 

70  °C 

temperature 

Recharge  procedure 

Charge  at  0.9  A  (CJl)  constant- 
current  rate  to  a  voltage  of  4. 1  V; 
continue  to  apply  a  constant- 
voltage  of  4.1  V  for  2.5  h  total 
recharge  time.  All  recharging  is 
to  begin  at  25  =t  3  °C 
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2.3.  Electrical  performance  tests 

Characterization  tests  were  performed  on  all  the  cells 
following  a  pre-test  readiness  review.  The  characterization 
tests  included  a  Cfl  static  capacity  test  (some  of  the  cells 
also  underwent  a  Cf 25  discharge  test);  low-  and  medium- 
current  hybrid  power  pulse  characterization  (L-HPPC  and 
M-HPPC,  respectively)  tests  for  18  s  at  discharge  currents  of 
2.7  and  7.2  A,  respectively;  and  a  7-day  self-discharge  test  at 
3.660  V,  which  corresponds  to  50%  SOC.  Thermal  perfor¬ 
mance  tests  consisting  of  the  static  capacity  and  L-HPPC 
tests  were  performed  on  four  cells  at  ambient  temperatures 
of  +5  and  +40  °C.  Finally,  reference  performance  tests 
(RPTs)  were  conducted  on  all  cells  prior  too  beginning 
the  life  testing.  The  RPTs  consisted  of  a  single  Ci/1  con¬ 
stant-current  discharge,  one  M-HPPC  test  and  ac  impedance 
measurements  at  1  kHz  at  100  and  0%  SOC.  Cells  were  then 
designated  for  either  the  calendar-  or  cycle-life  at  a  given 
condition.  The  cycle-  and  calendar-life  tests  where  then 
conducted  for  a  4-week  period  for  cells  tested  at  40,  50 
or  60  °C,  and  for  a  2-week  period  for  those  cells  tested  at 
70  °C.  The  RPTs  were  repeated  every  4  weeks  for  the  cells 
tested  at  40,  50  and  60  °C,  and  every  2  weeks  for  the  cells 
tested  at  70  °C.  End-of-testing  (EOT)  was  reached  when  a 
cell  was  unable  to  perform  the  M-HPPC  test  at  60%  DOD  or 
falling  below  the  3.0  V  minimum  cell  voltage. 

This  paper  will  deal  with  the  calendar-  and  cycle-life 
testing  of  these  cells  using  special  tests  developed  by  the 
ATD  program  [2].  The  special  calendar-life  test  is  theore¬ 
tically  charge-neutral,  so  it  will  not  perturb  the  SOC  of  the 
cell  any  more  than  absolutely  necessary.  The  magnitude  and 
duration  of  the  discharge  and  regen  tests  are  relatively 
modest,  compared  with  the  corresponding  M-HPPC  test 
profile,  so  they  would  have  a  minimal  affect  on  the  thermal 
condition  of  the  cell  under  test.  The  magnitude  and  duration 
of  the  discharge  pulse  for  this  test,  i.e.  3.6  A  for  9  s  was  set  to 
one-half  the  corresponding  values  used  for  the  M-HPPC  test. 
The  test  incorporated  a  somewhat  longer-than-normal  rest 
period  after  the  discharge  pulse  to  allow  additional  time  for 
voltage  recovery  before  the  regen  pulse.  This  is  a  new 
calendar-life  test  designed  to  obtain  additional  resistance 
data  at  regular  intervals  without  unduly  cycling  the  cells. 
The  idea  was  to  apply  a  single  pulse  profile  once  per  day 
from  which  the  discharge  and  regen  resistances  could  be 
calculated.  The  calendar-life  test  profile  is  shown  in  Table  1 . 


Table  1 

ATD  program  special  calendar-life  test  pulse  profile 


Step 
time  (s) 

Cumulative 
time  (s) 

Current 

(A) 

Charge 
(A  s) 

Cumulative 
charge  (A  s) 

9 

9 

3.6 

32.40 

32.40 

60 

69 

0.00 

0.00 

32.40 

2 

71 

-3.6 

-7.20 

25.20 

2 

73 

0.00 

0.00 

25.20 

47 

120 

-0.54 

25.38 

0.18 

Note  that  positive  values  for  the  current  correspond  to  a 
constant-current  discharge  and  negative  values  correspond 
to  a  constant-current  charge,  i.e.  regen  pulse. 

Each  cell  tested  using  the  calendar-life  test  was  assigned 
a  temperature  (40,  50,  60  or  70  °C)  and  a  target  SOC  (either 
60  or  80%  SOC  in  this  study).  For  a  new  cell,  the  cell  voltage 
was  3.742  Vat  60%  SOC  and  was  3.918  Vat  80%  SOC.  The 
value  of  the  voltage  at  a  given  SOC  was  determined  from  a 
calibration  table  provided  by  ANL  that  showed  the  voltage 
at  a  given  SOC,  as  found  by  conducting  Cf25  discharge 
tests  on  a  number  of  new  test  cells.  The  discharge  and 
regen  resistances  were  calculated  using  R  =  AV/AI,  i.e. 
the  change  of  the  voltage  of  the  cell  at  the  beginning  of 
the  discharge  (or  charge)  to  the  end  of  the  discharge  (or 
charge)  divided  by  the  change  in  the  current  during  the 
discharge  (or  charge).  For  the  calendar-life  tests  (Table  1), 
the  discharge  or  regen  currents  were  held  constant,  but  the 
voltage  did  change  during  the  course  of  the  discharge  or 
recharge.  The  test  was  conducted  once  per  day  for  a  4- week 
period  for  the  40,  50  and  60  °C  tests,  and  for  a  2- week  period 
for  the  70  °C  test. 

The  special  ATD  program  cycle-life  test  profiles  (con¬ 
ducted  at  60%  SOC)  for  the  A3,  A6  and  A9%  SOC  tests  are 
given  in  Tables  2-4,  respectively.  Positive  values  of  the 
current  correspond  to  a  discharge  current,  while  negative 
values  correspond  to  a  charge  current,  i.e.  regen.  These 
profiles  are  charge  neutral,  as  shown.  The  test  profiles  were 
conducted  once  the  cell  had  reached  the  test  temperature 

Table  2 


ATD  program  special  cycle-life  A3%  SOC  pulse  test  profile 


Step 
time  (s) 

Cumulative 
time  (s) 

Current 

(A) 

Charge 
(A  s) 

Cumulative 
charge  (A  s) 

14 

14 

7.20 

100.80 

100.80 

10 

24 

0.00 

0.00 

100.80 

2 

26 

-6.48 

-12.96 

87.84 

2 

28 

0.00 

0.00 

87.84 

32 

60 

2.745 

-87.84 

0.00 

20 

80 

Table  3 

ATD  program  special  cycle-life  A6%  SOC  pulse  test  profile 

Step 
time  (s) 

Cumulative 
time  (s) 

Current 

(A) 

Charge 
(A  s) 

Cumulative 
charge  (A  s) 

14 

14 

7.20 

100.80 

100.80 

10 

24 

0.00 

0.00 

100.80 

2 

26 

-6.48 

-12.96 

87.84 

2 

28 

0.00 

0.00 

87.84 

14 

42 

7.2 

100.80 

188.64 

10 

52 

0.00 

0.00 

188.64 

2 

54 

-6.48 

-12.96 

175.68 

2 

56 

0.00 

0.00 

175.68 

64 

120 

2.745 

-175.68 

0.00 

40 

160 
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Table  4 

ATD  program  special  cycle-life  A9%  SOC  pulse  test  profile 


Step 
time  (s) 

Cumulative 
time  (s) 

Current 

(A) 

Charge 
(A  s) 

Cumulative 
charge  (A  s) 

14 

14 

7.20 

100.80 

100.80 

10 

24 

0.00 

0.00 

100.80 

2 

26 

-6.48 

-12.96 

87.84 

2 

28 

0.00 

0.00 

87.84 

14 

42 

7.20 

100.80 

188.64 

10 

52 

0.00 

0.00 

188.64 

2 

54 

-6.48 

-12.96 

175.68 

2 

56 

0.00 

0.00 

175.68 

14 

70 

7.20 

100.80 

276.48 

10 

80 

0.00 

0.00 

276.48 

2 

82 

-6.48 

-12.96 

263.52 

2 

84 

0.00 

0.00 

263.52 

96 

180 

2.745 

-263.52 

0.00 

60 

240 

(either  40,  50,  60  or  70  °C).  Each  cell  undergoing  cycle- 
life  testing  was  initially  tested  at  the  target  temperature 
and  60%  SOC  for  100  iterations  followed  by  a  1  h  rest  to 
verify  stable  SOC  cycling.  If  the  SOC  cycling  was  stable, 
then  the  cycle-life  test  was  resumed  on  a  continuous  basis 
for  the  remaining  test  period.  Designated  cells  underwent 
the  cycle-life  testing  for  a  4-week  period  for  the  cells 
tested  at  40,  50  and  60  °C,  and  for  a  2-week  interval  for 
different  cells  tested  at  70  °C.  C/1  and  M-HPPC  reference 
tests  (at  25  °C)  were  performed  before  and  after  the  4-  or  2- 
week  test  interval. 


3.  Results  and  discussion 

3.1.  Calendar-life  tests  at  80%  SOC 

3.1.1.  Discharge  resistance 

For  the  calendar-life  tests  conducted  at  80%  SOC  (open- 
circuit  voltage  set  at  3.918  V),  cells  were  tested  at  both  the 
INEEL  and  SNL.  The  data  discussed  here  combined  all  of 
the  available  discharge  and  regen  resistances  as  a  function  of 
test  time  as  measured  by  each  of  the  two  laboratories.  Note, 
however,  that  some  of  the  cells  tested  under  a  given  tem¬ 
perature  condition  failed  prior  to,  or  during,  the  calendar-life 
test.  When  this  occurred,  the  data  for  that  cell  was  not 
included  in  the  averaged  data  presented  below.  The  notation 
convention  used  in  this  paper  to  identify  the  test  condition  is 
(xxCyy),  where  xv  gives  the  percentage  SOC  (either  60  or 
80%)  of  the  cell  during  the  test,  C  signifies  that  it  is  a 
calendar-life  test  and  yy  specifies  the  test  temperature  (either 
40,  50,  60  or  70  °C).  Thus,  (80C40)  specifies  a  cell  tested  at 
80%  SOC,  using  the  calendar-life  test  and  tested  at  40  °C. 
Fig.  1  shows  the  average  of  the  INEEL  and  SNL  (up  to  six 
cells)  discharge  resistances  at  40,  50,  60  and  70  °C  as  a 
function  of  the  square-root  of  the  time  at  test  temperature  in 
days.  The  data  shown  is  from  the  first  4- week  test  period  for 
the  40,  50  and  60  °C  tests,  and  for  the  first  2- week  test  period 
for  the  70  °C  test.  In  order  to  formulate  a  model  of  the 
behavior  of  the  resistance  data,  a  number  of  different  time 
dependencies  of  the  discharge  and  regen  resistances  were 
tried,  but  the  square-root  of  test  time  was  found  to  correlate 
the  data  the  best.  The  functional  forms  examined  were: 
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Fig.  1.  Average  of  all  INEEL  and  SNL  data  from  calendar-life  discharge  resistance  tests  at  (80C40,  80C50,  80C60  and  80C70)  for  ATD  Gen  1  cells  plotted  as 
a  function  of  the  square-root  of  test  time.  The  model  fits  to  the  relation:  R(t ,  T)  =  A(T)t +  B(T)  and  the  correlation  coefficient  ( R 2)  are  given  in  the  figure 
for  each  test  temperature. 
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R  =  At  +  B;  R  =  At2  +  Bt  +  C;  R  =  A[exp(Zfr)];  7?  = 
A[ln(7fr)]  +  C;  R  =  At3/2  +  7?;  R  =  A[^],  where  ^  is  a  fitting 
parameter  for  the  test  time  raised  to  a  power.  Higher  order 
polynomial  fits  to  the  data  would,  of  course,  give  better  fits 
to  the  data,  but  the  physical  significance  would  be  very  hard 
to  explain  in  terms  of  a  physical/chemical  model  of  the 
process(es)  (refer  to  [5]  for  further  details).  Also  shown  in 
the  Fig.  1  are  the  linear  regression  fits  to  the  function 
7?dis  =  At1/2  +  B ,  where  A  and  B  are  the  fitting  parameters. 
The  goodness-of-fit  parameter  is  given  by  the  value  R  , 
which  is  often  referred  too  as  the  correlation  coefficient. 
R2  =  1  would  be  the  best  correlation  of  the  data  to  the  fitting 
function.  In  general,  the  discharge  resistances  decrease  as 
the  test  temperature  is  increased.  This  is  true  for  all  of  the 
data,  except  at  70  °C,  where  the  resistance  is  slightly  greater 
than  that  measured  at  60  °C.  One  would  expect  that  the 
resistance  at  this  temperature  would  follow  the  general  trend 
and  be  lower  than  at  the  other  temperatures.  The  unusual 
70  °C  data  behavior  probably  indicates  that  an  additional 
mechanism  (or  mechanisms)  responsible  for  cell  degrada¬ 
tion  is  (are)  occurring  at  this  temperature. 

A  functional  form  for  the  description  of  the  time,  tem¬ 
perature  and  SOC  dependence  of  the  discharge  and  regen 
resistance  was,  therefore,  assumed  to  be  as  follows: 

R(t ,  r,  soc)  =  A(r,  soc)F(t)  +  b(t ,  soc)  (i) 

where  t  is  the  time  at  test  temperature  (in  days),  T  the  test 
temperature  and  SOC  the  SOC  of  the  battery  at  the  start  of 
the  calendar-life  test.  A(T,  SOC)  and  B(T,  SOC)  are  assumed 
to  be  functions  of  the  temperature  and  the  SOC  only.  The 
function  F(t)  is  assumed  to  be  only  a  function  of  the  time  at 
test  temperature;  the  best  fit  of  the  resistance  data  being 
F(t)  =  t1/2.  The  following  results  will  concern  the  verifica¬ 
tion  of  this  relationship  and  to  find  the  functional  forms  for 
A(T,  SOC)  and  B(T,  SOC).  Once  these  functions  have  been 
determined,  then  a  physical/chemical  basis  for  the  func¬ 
tional  forms  is  attempted  using  the  fits  to  the  discharge  and 
regen  resistances  as  a  guide.  The  reason  for  plotting  the  data 
as  various  functions  of  test  time  was  to  try  and  determine  not 
only  the  time  dependence  of  the  resistance  increase,  but  also 
if  the  a  functional  form  for  the  time  dependence  can  be 
found,  to  ascertain  a  physical/chemical  process  that  would 
account  for  the  discharge  and  regen  resistance  increase  with 
time.  This  information  could  then  be  used  to  understand  the 
process(es)  responsible  for  the  cell  degradation  and  in  turn 
suggest  possible  changes  in  the  construction  of  the  cells.  An 
additional  aspect  of  the  determination  of  the  functional  form 
of  the  time-dependent  degradation  would  be  to  predict  the 
calendar-  and  cycle-life  of  the  cells  at  various  test  tempera¬ 
tures. 

What  are  some  of  the  mechanisms  responsible  for  the 
resistance  and  increase  in  the  resistance  of  a  lithium-ion 
battery?  Zhang  et  al.  [8]  and  Wang  et  al.  [9]  have  discussed 
the  possible  mechanisms.  The  total  resistance  of  the  carbon 
anode  and  the  metal  oxide  cathode  is  the  sum  of  the 
following  resistances:  (a)  electrolyte  solution;  (b)  surface 


layer;  (c)  anode  and  cathode  particle  to  particle  contact;  (d) 
anode  and  cathode  particle  contact  to  the  current  collector; 
and  (e)  charge  transfer.  The  interfacial  impedance  at  the 
discharged  state  has  been  found  to  be  larger  when  compared 
with  the  charged  state  for  both  the  carbon  and  metal  oxide 
electrodes.  Experimental  results  using  electrochemical 
impedance  spectroscopy  (EIS)  show  that  the  impedance 
of  lithium-ion  cells,  at  least  with  LiCo02  cathodes,  is 
dominated  by  the  positive  electrode,  i.e.  the  cathode.  The 
total  cell  impedance  was  found  to  increases  with  a  decrease 
in  the  SOC.  Upon  consideration  of  the  multitude  of  possible 
mechanisms  that  can  lead  to  resistance  increases  as  the  cell 
ages,  the  fact  that  a  thin  film,  often  referred  too  as  the  solid 
electrolyte  interface  (SEI)  layer  arising  from  the  decom¬ 
position  of  the  electrolyte  and  salt  is  a  likely  candidate. 
When  the  electrodes  are  in  the  charged  state,  there  will  be  a 
large  portion  of  Ni4+  and  Co4+  cations  present  in  the 
cathode.  These  ions  have  a  strong  oxidizing  power  and 
can  react  with  the  electrolyte  and  salt  at  the  cathode/elec- 
trolyte  interface.  This  reaction  can  cause  decomposition  of 
the  electrolyte  and  salt  to  form  a  SEI  layer  on  the  cathode. 
After  extended  cycling,  the  LiNi1_^Cox02  electrode  used  in 
the  ATD  program  Gen  1  cells  should  be  heavily  passivated, 
resulting  in  a  large  resistance  at  the  electrolyte/cathode 
interface.  Due  to  this  increase  in  resistance,  the  reaction 
rate  will  be  lower  for  both  lithium-ion  insertion  (intercala¬ 
tion)  and  removal  (de-intercalation). 

The  earliest  reference  to  the  SEI  layer  that  the  authors  are 
aware  of  is  that  by  Thomas  et  al.  [10],  who  state  that  the 
polymeric  surface  layer,  i.e.  the  SEI  layer  must  be  in  a 
dynamic  state  that  depends  on  cell  temperature,  SOC  and  the 
extent  of  aging  of  the  cell.  The  resistances  directly  relate  to 
the  thickness  of  the  surface  layer.  The  SEI  layer  that  forms 
on  carbonaceous  electrode  materials  consists  of  many  dif¬ 
ferent  materials,  including  LiF,  Li2C03,  LiCO-R,  Li20, 
lithium  alkoxides  (Li-O-R,  where  R  is  a  hydrocarbon), 
non-conductive  polymers  and  a  number  of  other  possible 
chemical  compounds  composed  of  electrolyte  and  salt 
decomposition  products.  The  formation  of  the  SEI  layer 
mainly  occurs  on  the  anode  during  the  initial  formation 
(charging)  cycle  of  the  battery.  The  implication  of  the  SEI 
layer  on  the  carbon  electrode  is  that  it  will  cause  a  voltage 
drop  across  the  layer.  This  will  in  turn  modify  the  structure 
of  the  double  layer  at  the  carbon  electrode/electrolyte  inter¬ 
face,  which  generally  increases  the  charge  transfer  resis¬ 
tance  at  this  interface.  A  SEI  layer  also  forms  on  the  cathode 
material  that  also  results  in  an  increase  in  the  charge  transfer 
at  the  electrolyte/cathode  interface.  Cycling  will  also  cause 
capacity  loss  due  to  damage  and  disorder  in  the  metal  oxide 
cathode  particles.  Cycling  can  also  induce  severe  strain,  high 
defect  densities  and  occasional  fracture  of  the  particles  [8]. 
Severely  strained  cathode  particles  exhibit  cation  disorder. 
These  processes  lead  to  changes  in  the  thermodynamic 
properties  and  contact  resistance  of  the  metal  oxide  parti¬ 
cles.  The  accumulation  of  strain  in  the  electrode  particles 
may  cause  partial  shedding  of  the  electrode  material  from  its 
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current  collector.  A  portion  of  the  lithium-ions  in  the  cathode 
can  also  become  inactive  due  to  cation  disorder.  However, 
the  main  loss  in  the  cathode,  for  example,  LiCo02  is  caused 
by  the  change  in  resistance  at  the  surface  of  the  particles. 

Arora  et  al.  [1 1]  have  also  discussed  some  of  the  processes 
known  to  result  in  capacity  fade  in  lithium-ion  cells.  These 
are  lithium  deposition  on  the  anode  (over-charge  condition), 
electrolyte  decomposition,  anode  and/or  cathode  active 
material  dissolution,  phase  changes  in  the  anode  and  cathode 
materials  and  passive  film  formation  over  the  electrode  and 
current  collector  surfaces  (SEI  layer  formation).  They  found 
that  the  ionic  conductivity  of  the  electrolyte  does  not  con¬ 
tribute  significantly  to  the  cell’s  overall  conductivity.  This 
was  substantiated  by  the  work  of  Narayanan  et  al.  [12]  who 
state  that  the  process  of  lithium-ion  diffusion  in  the  anode 
and  cathode  lattice  is  considerably  slower  that  that  in  the 
electrolyte.  Therefore,  the  lithium-ion  diffusion  in  the  elec¬ 
trode  materials  would  be  one  of  the  rate-limiting  steps.  Thus, 
the  processes  affecting  the  cell  impedance  directly  relate  to 
the  electrode  materials  and  their  interactions  with  the  elec¬ 
trolyte,  i.e.  the  SEI  layer.  Ozawa  [13]  and  Megahed  and 
Scrosati  [14]  have  also  discussed  and  confirmed  these 
processes.  Nagasubramanian  et  al.  [15]  and  Baker  et  al. 
[16],  using  EIS  methods  find  that  the  impedance  increase 
with  aging  is  mostly  due  to  the  cathode.  They  also  found  that 
the  interfacial  impedance  increases  as  the  SOC  of  the  cell 
decreases.  From  their  measurements,  they  find  that  the  cell 
impedance  comes  mostly  from  the  cathode/electrolyte  inter¬ 
face,  not  from  the  anode/electrolyte  interface.  Guyomard 
and  Tarascon  [17]  conclude  that  oxidation  of  the  electrolyte 
is  the  main  failure  mechanism  for  lithium-ion  batteries. 

The  extensive  work  of  Auerbach  and  co-workers  [18,19] 
has  also  given  an  overview  of  the  processes  occurring  in  a 
lithium-ion  cell.  They  state  that  in  parallel  to  the  flux  of 
lithium-ions  to  and  into  the  electrodes,  there  is  a  flux  of 
electrons  from  the  current  collector  to  the  anode  or  cathode 
materials,  which  balances  the  charge.  This  electron  flux  also 
has  to  overcome  the  resistance  that  exists  among  the  elec¬ 
trode  particles,  all  of  which  are  partially  covered  by  elec¬ 
tronically  insulating  surface  films,  i.e.  a  SEI  layer.  Lithium 
intercalation  into  the  graphite  anode  or  the  metal  oxide 
cathode  is  a  serial  multi-step  process  in  which  lithium-ions 
have  to  first  migrate  through  the  electrolyte  and  then  through 
the  surface  films  that  covers  the  electrodes.  After  this 
migration,  insertion  into  the  electrode  material  is  accom¬ 
panied  by  a  charge  transfer  at  the  film/electrode  material 
interface.  This  is  then  followed  by  solid  state  diffusion  of 
lithium  into  the  bulk  of  the  electrode  material.  Finally, 
lithium  accumulates  within  crystallographic  sites  in  the 
bulk  electrode  material  via  phase  transition(s)  between 
the  various  intercalation  stages.  The  intercalation  stages, 
particularly  for  the  metal  oxide,  depend  on  the  crystalline 
structure  of  the  electrode  material.  The  process  of  charge 
transfer  resistance  can  be  related  to  three  different  processes: 

(1)  Li-ion  transfer  at  the  solution- surface  film  interface; 

(2)  Li-ion  transfer  at  the  surface  film/electrode  interface; 


and  (3)  inter-particle  electron  transfer  between  the  particles 
constituting  the  electrode  material.  They  also  state  that  the 
increased  resistance  observed  upon  cycling  the  battery 
mostly  reflects  changes  in  the  surface  structure  of  the 
electrodes.  After  prolonged  cycling,  there  are  phenomena 
such  as  expansion  and  contraction  of  the  electrode  material’s 
volume,  which  leads  to  local  breakdown  of  the  passivation 
layer  of  the  electrode  (on  a  microscopic  level).  This  allows 
continuous  reduction  and  oxidation  of  the  electrolyte/solute 
species.  While  the  process  occurs  on  a  very  small  scale,  it 
thickens  the  surface  films  and  consequently  the  electrode’s 
impedance  increases,  particularly  in  the  time  constants  that 
relate  to  lithium-ion  migration  through  the  surface  films, 
whose  increasing  thickness  upon  cycling  makes  them  more 
resistive.  The  electrolyte  composition  has  a  great  impact  on 
the  surface  films  and,  depending  on  its  composition,  the 
surface  films  may  be  the  dominant  factor  that  determines  the 
impedance  of  the  electrodes.  However,  this  behavior  may 
not  be  stable,  i.e.  the  electrode’s  impedance,  especially  in  the 
features  that  relate  to  the  surface  films,  increases  upon 
storage,  and  may  also  change  as  a  result  of  thermal  cycling 
and  the  charging  and  discharging  of  the  battery. 

In  summary,  the  overall  insertion  process  of  lithium  into 
the  battery  electrodes  is  quite  complicated.  It  includes 
diffusion  of  lithium-ions  in  the  solution  phase,  their  migra¬ 
tion  through  the  surface  films  (SEI  layer)  covering  the 
electrode  particles  (which  are  ionically  conducting  but 
electrically  insulating  to  electron  flow),  solid  state  diffusion, 
accumulation/consumption  of  lithium  in  the  bulk  (accom¬ 
panied  by  a  flux  of  electrons  that  counterbalance  the  charge), 
and  finally  phase  transition(s)  among  the  crystalline  struc¬ 
tures  of  the  electrode  materials.  Thus,  a  lithium-ion  battery 
is  a  very  dynamic  system  whose  electrical  performance 
depends  on  its  construction,  the  materials  used  in  its  assem¬ 
bly,  the  rates  of  charge  and  discharge  during  use,  the  SOC, 
and  its  temperature.  One  physical/chemical  process  that 
stands  out  as  a  candidate  for  having  the  greatest  impact 
on  the  impedance  of  the  cell  is  the  SEI  layer,  its  growth, 
composition,  structure  and  thickness.  These  attributes  will 
also  change  as  the  cell  ages. 

Looking  at  possible  analogous  processes  that  grow  thin 
films  upon  a  solid  surface,  one  can  consider  the  oxidation  of 
metals  [19-21].  Upon  examination  of  the  various  reaction 
rates  and  corresponding  rate  equations  for  the  oxidation  of 
metals,  it  is  found  that  they  are  functions  of  a  number  of 
factors,  such  as  material  composition,  temperature,  oxygen 
pressure,  elapsed  time  of  reaction,  surface  preparation  and 
pre-treatment  of  the  metal.  Although  rate  equations  alone 
are  insufficient  for  interpretations  of  oxidation  mechanisms, 
these  equations  may  be  used  to  classify  the  oxidation  of 
metals  and  may  as  such  often  limit  the  interpretation  to  a 
class  of  alternative  mechanisms.  The  rate  equations  most 
commonly  encountered  may  be  classified  as  logarithmic, 
parabolic  and  linear.  They  represent  only  limiting  and  ideal 
cases.  Deviations  from  these  rate  equations  and  intermediate 
rate  equations  are  also  often  encountered.  In  many  instances, 
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it  may  be  difficult  to  fit  rate  data  to  any  simple  rate  equation 
or  combination  of  rate  equations.  In  the  following  discus¬ 
sion,  an  analogy  is  made  between  a  process  occurring  at  the 
various  surfaces  present  in  a  lithium-ion  battery,  the  exact 
nature  having  been,  as  yet,  not  definitively  determined  and 
the  growth  of  an  oxide  film  on  a  metal  surface  [19-22].  Of 
the  various  oxidation  processes  that  have  been  observed,  the 
parabolic  rate  process  has  been  identified  and  has  the  func¬ 
tional  form: 

X  oc  f1/2  (2) 


where  the  thickness  of  the  thin  film  (X)  is  proportional  to  the 
square-root  of  the  time  that  the  film  growth  has  occurred. 
Blomgren  [23]  has  used  this  argument  to  account  for  the 
capacity  fade  of  the  lithium-ion  batteries  that  he  has  studied. 
As  a  rule,  parabolic  oxidation  signifies  that  a  thermal 
diffusion  process  is  rate  determining  [19,20].  Thermal  diffu¬ 
sion  processes  generally  have  a  temperature  dependence 
given  by  an  Arrhenius-like  process  where  the  diffusion  is 
given  by  [24-27]: 


exp 


where  D  is  the  diffusion  constant  in  cm  /s,  D0  the  limiting 
diffusion  constant  at  very  high  temperature,  E  the  activation 
energy  associated  with  the  diffusion  process,  R  the  gas 
constant  and  T  the  temperature  (K).  Such  a  process  may 
include  a  uniform  diffusion  of  one  or  both  of  the  reactants 
through  a  growing  scale  or  a  uniform  diffusion  of  gas  into 
the  metal. 

The  analogy  to  the  case  of  the  lithium-ion  battery  is  that 
there  would  be  a  film  growing  on  the  surface  of  the  anode 
and/or  cathode  materials  over  a  period  of  time  that  would  be 
temperature-  and  possibly  SOC-dependent.  The  nature  of 
the  thin  film  would  also  depend  on  the  electrolyte  and  the 
composition  of  the  electrodes.  The  thickness  of  this  thin  film 
could  give  rise  to  an  increase  in  the  resistance  of  the  cell  as 
the  rate  of  lithium  migration  into/out  of  the  anode  and/or 
cathode  materials  would  be  impeded  by  the  thin  film.  The 
thicker  the  thin  film,  the  lower  the  mobility  of  the  lithium- 
ions  and  thus  the  higher  the  resistance. 

The  best  fit  to  the  time  dependence  of  the  resistance,  as 
previously  discussed,  is  the  square-root  of  the  time  at  test 
temperature  as  shown  in  Fig.  1.  It  may  well  be  that  as  the  cell 
ages  the  SEI  layer  grows  in  thickness,  leading  to  an  increase 
in  the  resistance,  due  to  a  decrease  in  the  migration  rate  of 
the  lithium-ions  into/out  of  the  anode  and/or  the  cathode. 
The  stresses  experienced  by  the  cathode  particles  during 
change  and  discharge,  and  during  temperature  variation, 
could  lead  to  fracturing  of  the  particles,  as  mentioned.  This 
would  effectively  expose  new  surfaces  on  which  a  SEI  layer 
would  grow,  thus  increasing  the  anode  and/or  cathode 
resistance.  This  resistance  would  be  observed  as  an  increase 
of  the  discharge  and  charge  resistances  measured  during  the 
calendar-life  test.  If  the  increase  in  the  discharge  and  regen 
resistance  is  proportional  to  the  square-root  of  time  at  test 


temperature,  then  the  resistances  can  be  expressed  by  a 
function  having  the  form  (from  Eq.  (1)): 

R(t,  T,  SOC)  =  A(T,  SOC )tl/2  +  B(T,  SOC)  (4) 


where  the  discharge  and  regen  resistance  ( R(t ,  T \  SOC))  is  a 
function  of  time  (t),  test  temperature  (7)  and  SOC.  That  there 
is  a  dependence  on  the  SOC  is  verified  by  comparing  the 
discharge  and  regen  resistance  when  the  calendar-life  test  is 
conducted  at  80  or  60%  SOC  (discussed  later  in  this  paper). 
The  functions  A(T,  SOC)  and  B(T ]  SOC)  are  assumed  to  be 
functions  of  the  test  temperature  and  SOC.  To  determine  the 
temperature  dependence  of  the  functions  A(T,  SOC)  and  B(T, 
SOC),  one  can  plot  the  fitting  coefficients  determined  from 
the  fits  shown  in,  for  example,  Fig.  1  as  various  function  of  the 
test  temperature.  Shown  in  Fig.  2  is  a  plot  of  the  function  A(T, 
SOC)  for  the  discharge  resistance  determined  from  the  fits 
shown  in  Fig.  1  as  a  function  of  the  inverse  of  the  test 
temperature  (K).  Other  functional  dependencies  were  exam¬ 
ined  [5],  but  this  relation  fit  the  A(T,  SOC)  parameter  the  best. 
In  Fig.  2,  the  A(T,  SOC)  parameter  is  plotted  as  an  inverse 
function  of  the  temperature  (K)  using  an  exponential  fit,  i.e. 
A(T,SOC)  =  a(SOC){exp[/?(SOC)/7]}  where  the  fitting 
parameters  a  and  b  are  shown  as  being  function  of  the 
SOC  of  the  battery.  This  type  of  function  occurs  in  a  number 
of  physical/chemical  processes,  such  as  chemical  kinetics  and 
in  diffusion.  A  functional  form  of  this  kind  is  generally 
referred  too  as  Arrhenius-like  behavior,  as  it  applies  to  a 
chemical  reaction  or  to  a  diffusion- type  of  process.  The  fit  of 
A(T,  SOC)  to  this  functional  form  of  the  temperature  is  quite 
good  for  40,  50  and  60  °C,  but  the  70  °C  value  does  not  fall  on 
the  fitting  curve.  The  fits  to  the  exponential  function  given  in 
Fig.  2  are  given  for  when  only  the  40,  50  and  60  °C  data  are 
used,  and  also  for  when  all  the  temperature  data  are  used,  i.e. 
the  40,  50,  60  and  70  °C  data.  The  values  of  the  fitting 
parameters  and  the  R '  value  are  given  in  the  figure.  Using, 
then,  the  fit  to  the  exponential  of  the  inverse  temperature,  the 
A  parameter  can  be  expressed  as: 


A(T,SOC)  =  a(SOC)<  exp 


'fc(SOC)' 


T 


where  a( SOC)  is  the  fitting  coefficient  having  the  units  of 
(Q/t172)  and  b( SOC)  has  the  units  of  temperature  (K).  The  b 
fitting  parameter  can  be  related  to  an  activation  energy 
(EdisA)  using  b  =  Edis^/R  where  R  is  the  gas  constant  equal 
to  8.314  J/(mol  K)  (or  1.987  cal/(mol  K)). 

Similarly,  the  discharge  B(T,  SOC)  parameters  from  the 
fits  to  the  data  shown  in  Fig.  1  are  plotted  as  a  function  of  the 
inverse  of  the  test  temperature  (K)  (Fig.  3).  Two  fits  are  again 
shown,  one  for  when  only  the  40,  50  and  60  °C  data  are  used 
and  one  for  when  all  of  the  temperatures  are  used.  Again,  the 
70  °C  value  does  not  correlate  well  with  the  other  tempera¬ 
tures.  The  most  reasonable  fit  to  the  discharge  B(T,  SOC) 
parameter  was  thus  found  to  be  given  by  the  expression: 


B(T,  SOC)  =  c( SOC) 


exp 


'd(SOC)' 


T 
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[1/Test  Temperature]  (1/K) 


O  Temperatures  40,50,60,70  °C  ♦  Temperatures  40,50,60  °C 
—  Exp  (Temperatures  40,50,60  °C)  -  -  -  Exp  (Temperatures  40,50,60,70  °C) 


Fig.  2.  A(T)  parameter  for  the  tests  (80C40,  80C50,  80C60  and  80C70)  from  data  fits  to  the  model  R(t ,  T)  =  A(7)A2  +  B(T)  is  shown  as  a  function  of  the 
inverse  test  temperature  (K).  Fitting  parameters  to  the  model  A (T)  =  a[&xp(b/T )]  are  given.  The  fits  to  when  all  the  test  temperatures  are  used  and  when  only 
the  40,  50  and  60  °C  data  are  used  are  shown. 


where  c(SOC)  is  a  coefficient  having  the  units  in  Q  and 
d( SOC)  is  the  fitting  parameter  having  the  units  of  tem¬ 
perature  (K).  As  before,  the  c  and  d  parameters  are  shown 
to  possibly  be  functions  of  the  SOC  of  the  battery.  The  b 
parameter  can  be  related  to  an  activation  energy  using 


b  =  EdisB/R,  where  R  is  again  the  gas  constant.  Other 
functions  for  the  temperature  dependence  of  A(T,  SOC) 
and  B(T ]  SOC)  were  examined,  but  the  exponential  of 
the  inverse  of  the  test  temperature  was  found  to  be  the 
best  [5]. 


[1/Test  Temperature]  (1/K) 


O  Temperatures  40,50,60,70  °C  •  Temperatures  40,50,60  °C) 

- Exp  (Temperatures  40,50,60  °C)  -  -  -  Exp  (Temperatures  40,50,60,70  °C) 


Fig.  3.  B(T)  parameter  for  the  tests  (80C40,  80C50,  80C60  and  80C70)  from  data  fits  to  the  model:  R(t ,  T)  =  A(7)A2  +  B(T)  is  shown  as  a  function  of  the 
inverse  test  temperature  (K).  Fitting  parameters  to  the  model  B(T)  =  c[exp(d/7)]  are  given.  The  fits  to  when  all  the  test  temperatures  are  used  and  when  only 
the  40,  50  and  60  °C  data  are  used  are  shown. 
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Table  5 

Values  of  Eact  at  80%  SOC  from  analysis  of  calendar-life  test  data  using  the  relationship:  R(t,  T)  =  A(T)tl/2  +  B(T),  where  A(T)  =  a[exp(EactA/RT)\  and 
B(T)  =  c[exp(EactB/RT)] 


Test  condition 

Activation  energy  (Ea cuA,  kJ/mol) 

Activation  energy  (Eact  B,  kJ/mol) 

Rdis  (80C40,  80C50,  80C60  and  80C70) 

1.58  (0.378) 

9.04  (2.16) 

Rdis  (80C40,  80C50  and  80C60) 

7.74  (1.85) 

12.2  (2.91) 

tfreg  (80C40,  80C50,  80C60  and  80C70) 

4.36  (1.04) 

9.20  (2.20) 

Rre g  (80C40,  80C50  and  80C60) 

8.03  (1.92) 

12.9  (3.08) 

Eact  are  activation  energies,  and  a  and  c  are 

pre-exponential  constants.  The  values  given  in  parenthesis 

are  in  kcal/mol. 

Table  6 

Values  of  the  pre-exponential  constants  a  and  c  at  80%  SOC  from  analysis  of  calendar-life  test  data  analysis  using  the  relationship:  R(t,  T)  =  A(T)tl/2  +  B(T), 
where  A(T)  =  a[exp(EactA/RT )]  and  B(T)  =  c[exp(EacUB/RT)] 

Test  condition 

Pre-exponential  constant  a  (x  10  4  Q  per  (day)172) 

Pre-exponential  constant  c  (x  10  4  Q) 

Rdis  (80C40,  80C50  and  80C60) 

1.29 

4.47 

RIeg  (80C40,  80C50  and  80C60) 

1.03 

2.68 

The  values  for  the  activation  energies  and  the  pre-expo¬ 
nential  constants  determined  from  the  fits  to  A(T,  SOC)  and 
B(T,  SOC)  for  both  the  discharge  and  regen  resistance  are 
given  in  Tables  5  and  6.  Note  that  the  parameters  a,  b ,  c  and  d 
can  be  SOC-dependent,  which  is  indeed  the  case,  as  will  be 
shown  when  the  calendar-life  60%  SOC  discharge  and  regen 
resistance  data  are  analyzed.  However,  it  should  be  restated 
that  since  the  calendar-life  tests  were  only  conducted  at  two 
SOC,  a  model  for  the  SOC  dependence  could  not  be 
determined,  as  only  two  data  points  were  available.  The 
final  relationship  for  the  discharge  (and  regen)  resistance  is, 


therefore,  the  following: 


R(t ,  r,SOC)  =  a(SOC)\  exp 


+  c(SOC)  <  exp 


X  soc)' 

T 

d(SOC) 


t 


1/2 


T 


Fig.  4  shows  the  model’s  predictions  using  Eq.  (7)  for  the  case 
when  all  of  the  temperature  values  (40,  50,  60  and  70  °C)  are 
used  in  determining  a ,  b ,  c  and  d  for  the  discharge  resistance. 
The  raw  data  values  are  the  same  as  shown  in  Fig.  1  and  are 
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Fig.  4.  Model  predictions  compared  to  calendar-life  discharge  resistance  data  for  ATD  Gen  1  (80C40,  80C50,  80C60  and  80C70)  tests  with  the  a,  b,  c  and  d 
parameters  determined  when  all  the  test  temperatures  were  used  to  fit  the  resistance  data. 
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Fig.  5.  Model  predictions  compared  to  calendar-life  discharge  resistance  data  for  ATD  Gen  1  (80C40,  80C50,  80C60  and  80C70)  tests  with  the  a ,  b,  c  and  d 
parameters  determined  when  only  the  40,  50  and  60  °C  test  temperatures  were  used  to  fit  the  resistance  data.  The  model  prediction  for  a  test  temperature  of 
25  °C  is  also  shown. 


presented  with  the  model  predictions.  The  fits  are  not  bad,  but 
the  70  °C  data  are  not  well  correlated,  as  should  be  the  case,  as 
previously  discussed.  Fig.  5  shows  fits  to  the  discharge 
resistance  data  using  the  model  given  in  Eq.  (7)  when  the 
values  for  a ,  b ,  c  and  d  were  determined  using  only  the  40,  50 
and  60  °C  resistance  data.  This  model  gives  excellent  fits  to 
the  40,  50  and  60  °C  resistance  data.  The  model  using  these 
fitting  parameters  does  not  fit  the  70  °C,  as  would  be  expected. 
Also  included  in  this  figure  are  the  model  predictions  for  the 
discharge  resistance  as  a  function  of  time  at  test  temperature 
when  the  test  temperature  is  25  °C. 

3.1.2.  Regen  resistance 

This  section  treats  the  regen  resistance  for  the  calendar- 
life  testing  at  80%  SOC  of  the  ATD  Gen  1  cells.  The  same 
approach  was  used  as  was  used  for  the  treatment  of  the 
discharge  resistance  at  80%  SOC.  All  of  the  regen  resistance 
data  for  the  four  test  temperatures  for  all  the  INEEL  and 
SNL  tested  cells  is  given  in  Reference  [5]  and  will  not  be 
reproduced  here.  Fig.  6  shows  the  average  of  all  of  the 
available  regen  resistance  data  at  each  of  the  four  test 
temperatures  for  the  first  4- week  period  at  40,  50  and 
60  °C,  and  for  the  first  2- week  period  at  70  °C.  As  was 
the  case  for  the  discharge  resistance,  the  regen  resistance 
increased  non-linearly  with  increasing  time  at  test  tempera¬ 
ture.  The  regen  resistance  was  smaller  the  higher  the  test 
temperature,  as  was  also  the  case  with  the  discharge  resis¬ 
tance.  The  regen  resistance  was  also  smaller  than  the  dis¬ 
charge  resistance  at  each  of  the  respective  test  temperatures. 
The  70  °C  resistance,  as  shown  in  Fig.  6,  did  not  follow  the 
general  trend  of  the  other  temperatures  in  that  the  resistance 


was  the  same  or  higher  than  the  60  °C  temperature  resis¬ 
tance.  This  trend  for  the  70  °C  regen  data  is  the  same  as  was 
observed  for  the  70  °C  discharge  resistance  (Fig.  5).  This 
indicates  again  that  something  has  happened  to  the  cells  at 
this  test  temperature  compared  to  the  lower  test  tempera¬ 
tures.  The  detailed  nature  of  the  process(es)  responsible  for 
this  increase  in  the  resistance  is  (are)  not  presently  known. 

The  same  approach  was  used  to  analyze  the  regen  resis¬ 
tance  at  the  various  temperatures.  The  relations  as  given  in 
Eqs.  (l)-(6)  were  also  found  to  best  correlate  the  time  and 
temperature  dependence  of  the  regen  resistance,  but  with 
different  values  of  the  fitting  parameters  a ,  b ,  c  and  d.  The 
best  correlation  to  the  time  dependence  of  the  regen  resis¬ 
tance  was  obtained  using  a  square-root  of  the  time  at  test 
temperature,  as  was  found  to  be  the  case  for  the  discharge 
resistance.  The  fits  to  this  time  dependence,  the  fitting 
parameters,  and  the  goodness-of-fit  values  are  also  shown 
in  Fig.  6.  The  fits  are  quite  good  for  all  of  the  test  tempera¬ 
tures.  Thus,  for  the  regen  resistance,  the  functional  form  for 
the  time  dependence  is  that  given  by  Eq.  (4),  which  was  also 
used  for  the  discharge  resistance.  The  fitting  parameters  for 
each  temperature,  i.e.  the  slope  being  A(T,  SOC)  and  the 
intercept  being  B(T,  SOC),  were  then  used  to  determine 
the  temperature  dependence  of  A  and  B  as  before.  As  was  the 
case  for  the  discharge  resistance,  the  most  meaningful 
temperature  dependence  for  A(T, \  SOC)  and  B(T,  SOC) 
was  to  use  an  Arrhenius-like  function  to  determine  the 
temperature  dependence  of  A(T,  SOC)  and  B(T,  SOC),  as 
given  by  Eqs.  (5)  and  (6).  The  values  for  the  fitting  coeffi¬ 
cients  are  given  in  Tables  5  and  6.  The  best  overall  correla¬ 
tion  to  the  regen  resistance  was  when  fits  to  the  70  °C  data 
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Fig.  6.  Average  of  all  INEEL  and  SNL  data  from  calendar-life  regen  resistance  tests  at  (80C40,  80C50,  80C60  and  80C70)  for  ATD  Gen  1  cells  plotted  as  a 


function  of  the  square-root  of  test  time.  The  model  fits  to  the  relation:  R(t,  T) 
each  temperature. 

were  excluded  from  the  model.  Fig.  7  shown  the  predictions 
of  the  model  using  Eq.  (7)  and  the  parameter  values  given  in 
Tables  5  and  6  for  when  only  the  40,  50  and  60  °C  data  fits 
are  used.  Also  shown  in  Fig.  7,  is  the  prediction  of  the  model 
when  the  test  temperature  is  25  °C.  The  model  predictions 


A(T)t1/2  +  B(T)  and  the  correlation  coefficient  (. R 2)  are  given  in  the  figure  for 


shown  in  Fig.  7  for  the  regen  resistance  can  be  compared  to 
the  model  predictions  shown  in  Fig.  5  for  the  discharge 
resistance.  Overall,  the  discharge  resistance  is  greater  than 
the  regen  resistance  at  all  times  at  test  temperature  and  at  all 
test  temperatures. 
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Fig.  7.  Model  predictions  compared  to  calendar-life  regen  resistance  data  for  ATD  Gen  1  (80C40,  80C50,  80C60  and  80C70)  tests.  Model  predictions  are 
shown  with  the  a,  b,  c  and  d  parameters  determined  when  only  the  40,  50  and  60  °C  test  temperatures  were  used  to  fit  the  resistance  data.  The  model 
prediction  for  a  test  temperature  of  25  °C  is  also  shown. 
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Table  5  shows  the  activation  energies  for  both  the  dis¬ 
charge  and  regen  resistance  data  as  determined  using  Eq.  (7) 
with  the  square-root  of  test  time  at  temperature  time  depen¬ 
dence  and  the  exponential  of  the  inverse  temperature  for  the 
temperature  dependence.  Given  in  this  table  are  the  activa¬ 
tion  energies  for  A(T,  SOC)  and  B(T. \  SOC)  for  when  the 
70  °C  data  were  used  in  the  data  fits  and  for  when  only  the 
40,  50  and  60  °C  data  were  used  to  determine  the  para¬ 
meters.  The  fits  to  the  40,  50  and  60  °C  data  are  the  best  ones 
to  consider  based  on  the  observations  discussed.  Upon 
examination  of  the  activation  energies  given  in  Table  5, 
general  trends  for  the  activation  energies  are  that  the  dis¬ 
charge  and  regen  values  are  approximately  the  same  order  of 
magnitude  for  EactA  and  for  EactB.  The  activation  energies 
for  EactB  are  larger  than  for  EactA.  The  activation  energies 
for  Rdis  and  Rreg  are  approximately  the  same  within  experi¬ 
mental  error  for  both  Eac t  A  and  EactB.  There  is  presently  no 
explanation  for  the  order  of  magnitude  of  these  values  or 
why  the  EactB  values  are  greater  than  the  EactyA  values. 
Fellner  et  al.  [28]  have  discussed  measurements  on  a  battery 
when  the  cathode  was  a  mixture  of  FiNi02  and  FiCo02,  the 
electrolyte  was  a  mixture  of  PC,  EC  and  DMC  with  LiPF6, 
and  the  anode  was  a  mixture  of  two  different  carbons.  They 
observed  that  the  interfacial  resistance  for  a  cell  was  much 
higher  at  lower  temperatures  and  that  the  interfacial  resis¬ 
tance  appeared  to  grow  linearly  with  test  cycle  number.  The 
interfacial  resistance  had  an  activation  energy  of  63.5  kJ/mol. 
This  value  for  the  activation  energy  was  attributed  to  an 
activated  rate  process  such  as  solid  state  diffusion  and/or 
electrochemical  kinetics.  The  present  authors  have  been 
unable  to  find  literature  values  for  lithium-ion  diffusion 
through  the  SEI  layer  on  the  anode  or  cathode,  or  its 
diffusion  in  the  anode  or  cathode  materials  for  the  composi¬ 
tions  used  in  the  ATD  Gen  1  cells. 

For  completeness,  the  values  of  the  pre-exponential  con¬ 
stants  obtained  from  the  fits  for  the  case  of  only  using  the  40, 
50  and  60  °C  data  in  the  fits  are  shown  in  Table  6.  The  values 
for  the  pre-exponential  parameters  a  and  c  are  greater  for  the 
discharge  resistance  than  for  the  regen  resistance.  At  present, 
there  is  no  explanation  for  these  values  based  on  a  physical/ 
chemical  model,  nor  have  literature  values  been  found  that 
may  provide  insight  into  the  values. 

3.2.  Calendar-life  tests  at  60%  SOC 

This  section  discusses  the  discharge  and  regen  resistance 
acquired  when  the  cells  were  at  60%  SOC.  The  test  pulse 
sequence  was  that  shown  in  Fig.  F  The  tests  were  conducted 
at  four  temperatures:  40,  50,  60  and  70  °C,  using  the  same 
testing  protocols  used  for  the  80%  SOC  tests  discussed 
previously. 

3.2.1.  Discharge  and  regen  resistance 

This  section  discusses  the  calendar-life  test  (see  Fig.  1  for 
the  discharge/charge  test  sequence)  when  the  cells  were  at 
60%  SOC  (open-circuit  voltage  of  3.747  V  for  new  cells). 


All  of  the  test  data  were  acquired  at  the  INEEF.  Initially, 
three  cells  were  tested  at  each  of  the  four  temperatures:  40, 
50,  60  and  70  °C.  The  same  methodology  was  used  to 
analyze  these  data  as  was  used  to  analyze  the  80%  SOC 
discharge  and  regen  resistance.  Further  details  of  the  pro¬ 
cedure  and  the  results  can  be  found  in  [5].  At  this  SOC,  it 
was  observed  that  there  was  an  onset  of  a  new  process  that 
affected  the  resistance  values  at  60  and  70  °C.  The  discharge 
and  regen  resistance  data  as  a  function  of  time  at  test 
temperature  are  shown  in  Figs.  8  and  9,  respectively.  The 
best  fits  to  the  data  using  the  model  given  in  Eq.  (7)  are  also 
shown.  Tables  7  and  8  present  the  fitting  parameters  deter¬ 
mined  from  the  model  for  the  60%  SOC  calendar-life  tests. 
As  can  be  observed  for  the  discharge  resistance,  the  model  to 
the  40  and  70  °C  are  quite  good.  However,  the  discharge 
resistance  at  60  °C  was  measured  to  be  greater  than  the 
50  °C  resistance.  The  model,  however,  does  bracket  the  data, 
but  predicts  that  the  50  °C  discharge  resistance  should  be 
greater  than  the  60  °C  discharge  resistance.  The  60  °C  regen 
resistance  data  is  smaller  than  the  50  °C  resistance  as  would 
be  expected,  but  still  appears  to  be  anomalous.  It  was  also 
assumed  that  the  60  °C  data  does  not  follow  the  general  trend 
of  the  other  temperatures  due  to  an  unusual  process  occur¬ 
ring  at  this  temperature.  There  was  no  great  improvement  [5] 
in  the  correlation  of  the  model  with  the  data  when  the  60  °C 
data  was  excluded  as  compared  to  when  all  the  test  tem¬ 
peratures  are  used  to  fit  A(T,  SOC)  and  B(J \  SOC)  in  the 
model.  This  effect  may  well  be  a  result  of  the  SOC  influence 
on  the  behavior  of  the  cell.  This  influence  may  cause  new 
processes  to  occur  at  different  temperatures,  depending  on 
the  SOC  of  the  cell.  As  noted  when  discussing  the  80%  SOC 
data,  it  may  also  be  the  case  that  the  70  °C  resistance  data  are 
also  unusual  due  to  the  same  or  similar  processes  responsible 
for  the  behavior  during  testing  at  this  SOC. 

The  values  for  the  activation  energies  derived  from  the  fits 
to  the  (60C40,  60C50,  60C60  and  60C70)  test  discharge  and 
regen  resistance  data  sets  are  shown  in  Table  7.  These  values 
compared  to  the  (80C40,  80C50  and  80C60)  test  values  (see 
Table  5)  are  of  the  same  order  of  magnitude  and  are  of 
comparable  value.  The  EactB  values  are  smaller  at  (60C40, 
60C50,  60C60  and  60C70)  than  they  are  at  (80C40,  80C50 
and  80C60).  This  difference  could  partially  account  for  the 
lower  discharge  and  regen  resistance  at  60%  SOC  than  at 
80%  SOC.  The  exponential  of  the  activation  energy  divided 
by  the  temperature  term,  however,  is  only  one  part  of  the 
total  resistance  formula  as  given  by  Eq.  (7).  One  also  has  to 
consider  the  pre-exponential  factors.  The  pre-exponential 
factors  as  determined  from  the  data  for  the  (60C40,  60C50, 
60C60  and  60C70)  tests  are  given  in  Table  8. 

3.3.  Cycle-life  discharge  and  regen  resistance  as  a 
function  of  test  time,  A  %  SOC  and  test  temperature 

The  cycle-life  test  at  60%  SOC  was  only  conducted  at  the 
INEEF.  The  following  data  combined  all  of  the  available 
discharge  and  regen  resistances.  Note,  however,  that  some  of 
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Fig.  8.  Model  predictions  compared  to  calendar-life  discharge  resistance  data  for  ATD  Gen  1  (60C40,  60C50,  60C60  and  60C70)  tests.  Model  used  the  a,  b,  c 
and  d  parameters  determined  when  all  the  test  temperatures  were  used  to  fit  the  resistance  data.  The  model  prediction  for  a  test  temperature  of  25  °C  is  also 
shown. 


the  cells  tested  under  a  given  temperature  and  A%  SOC 
condition  failed  prior  too,  or  during,  the  cycle-life  test. 
When  this  occurred,  the  data  for  that  cell  was  not  included 
in  the  averaged  data  presented  below.  The  notation  conven¬ 
tion  used  in  this  paper  to  identify  a  cycle-life  test  is  the  same 
as  used  for  the  calendar-life  except  that  Z  in  ( xxZyy )  is 


replaced  by  3,  6  or  9  and  signifies  that  it  is  a  cycle-life  test  at 
either  A3,  A6  or  A9%.  Thus,  (60340)  would  specify  a  cell 
tested  at  60%  SOC,  using  the  A3%  SOC  swing  cycle-life  test 
and  tested  at  40  °C. 

Similar  to  the  calendar-life  modeling  discussed  previou¬ 
sly,  Fig.  10  shows  the  average  discharge  (designated  Rdis) 


Fig.  9.  Model  predictions  compared  to  calendar-life  regen  resistance  data  for  ATD  Gen  1  (60C40,  60C50,  60C60  and  60C70)  tests.  Model  used  the  a,  b,  c  and 
d  parameters  determined  when  all  the  test  temperatures  were  used  to  fit  the  resistance  data.  The  model  prediction  for  a  test  temperature  of  25  °C  is  also 
shown. 
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Table  7 


Values  of  Eact  at  60%  SOC  from  analysis  of  calendar-life  test  data  using  the  relationship:  R(t,  T)  =  A(T)tl/2 
B(T)  =  c[exp(EactB/RT)] 

+  B(T),  where  A(T)  =  a[Qxp(EactA/RT)]  and 

Test  condition 

Activation  energy  (EactA,  kJ/mol) 

Activation  energy  ( Eact  B ,  kJ/mol) 

Rdis  (60C40,  60C50,  60C60  and  60C70) 

5.88  (1.41) 

8.11  (1.94) 

Rreg  (60C40,  60C50,  60C60  and  60C70) 

8.35  (2.00) 

8.59  (2.05) 

£aCt  are  activation  energies,  and  a  and  c 

are  pre-exponential  constants.  The  values  given  in  parenthesis  are  in 

kcal/mol. 

Table  8 

Values  of  the  pre-exponential  constants  a  and  c  at  60%  SOC  from  analysis  of  calendar-life  test  data  using  the  relationship:  R(t,  T)  =  A(T)t1/2  +  B(T),  where 
A(T)  =  <2[exp(£act)A/Er)]  and  B(T)  =  c[Qxp(EacUB/RT)] 

Test  condition 

1  /O 

Pre-exponential  constant  a  (Q  per  (day)  ) 

o 

Pre-exponential  constant  c  (x  10“  Q) 

Rdis  (60C40,  60C50,  60C60  and  60C70) 

1.88  X  1CT4 

2.03 

Rreg  (60C40,  60C50,  60C60  and  60C70) 

6.56  x  1(T5 

1.31 

and  regen  (designated  Rreg)  resistances  as  a  function  of  the 
square-root  of  time  at  test  temperature  at  60%  SOC,  A3% 
SOC  swing  and  at  test  temperatures  of  40,  50,  60  and  70  °C. 
Note  that  each  data  point  shown  represents  the  discharge  or 
regen  value  measured  after  100  cycle-life  cycles  had  been 
completed,  i.e.  between  the  data  points  shown,  100  test 
cycles  were  done.  Note,  as  pointed  out  later,  that  it  was 
assumed  that  all  tests  at  a  given  temperature  were  actually 
done  at  that  temperature.  There  was,  however,  some  varia¬ 
tion  in  the  test  temperatures  for  a  given  cell  at  a  nominally 
specified  target  test  temperature.  This  variation  in  the  test 
temperature  was  not  accounted  for  in  Fig.  10,  as  the  variation 
was  usually  only  several  °C.  This  temperature  variation  is, 


however,  addressed  in  the  detailed  treatments  of  the  data 
presented  below.  The  averaged  data  shown  in  Fig.  10  have 
also  been  corrected  for  the  time  it  took  the  cells  to  stabilize, 
which  is  why  all  data  start  at  approximately  time  =  0  h.  In 
general,  the  resistance  decreases  as  the  temperature  is 
increased.  This  is  true  for  all  data  except  that  at  70  °C, 
where  the  resistance  is  slightly  greater  than  that  measured  at 
60  °C.  One  would  expect  that  the  resistance  at  this  tempera¬ 
ture  would  follow  the  general  trend  and  be  lower  than  at  the 
other  temperatures.  The  unusual  behavior  of  the  70  °C  data 
probably  indicates  that  an  additional  mechanism  (or 
mechanisms)  responsible  for  cell  degradation  is  (are)  occur¬ 
ring  at  this  temperature.  This  was  also  found  to  be  the  case 
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Fig.  10.  Average  of  all  INEEL  cycle-life  discharge  and  regen  resistance  data  for  the  A3%  cycle-life  (60340,  60350,  60360  and  60370)  tests  as  a  function  of 
the  square-root  of  the  test  time. 
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for  the  calendar-life  tests  as  previously  discussed.  It  can  also 
be  observed  in  Fig.  10  that  the  discharge  resistance  increased 
faster  than  the  regen  resistance  at  each  test  temperature  as 
the  slopes  of  the  discharge  data  are  greater  than  the  regen 
data.  The  data  shown  also  do  not  correspond  to  a  period  of  4 
weeks  for  the  40,  50  and  60  °C  tests,  or  to  a  2-week  period 
for  the  70  °C  test.  Due  to  instrumentation  problems,  the  data 
collection  was  interrupted  at  the  longer  test  times.  However, 
the  tests  did  continue  for  the  full  4-  or  2-week  period.  With 
regard  to  the  plots  showing  the  resistance  as  a  function  of  the 
square-root  of  test  time,  other  time-dependent  fitting  func¬ 
tions  were  used  in  attempts  to  correlate  the  data  given  in  this 
figure  [6].  The  square-root  of  test  time  gave  the  best  fit  to  the 
resistance  as  was  also  found  to  be  the  case  for  the  calendar- 
life  tests  described  above. 

Fig.  1 1  displays  the  fits  of  the  average  Rdis  to  the  square- 
root  of  test  time  (h),  using  an  expanded  y-axis  scale,  for  the 
four  test  temperatures.  The  slopes  of  the  best  fit  to  this  time- 
dependent  function,  as  well  as  the  time-equal-to-zero  inter¬ 
cepts  are  given  in  the  figure.  The  goodness  of  fit  at  each  of 
the  test  temperatures,  R  ,  is  also  given.  The  slopes  and  the 
intercept  values  are  generally  seen  to  decrease  with  increas¬ 
ing  temperature.  The  R2  values  are  also  quite  good.  The 
exception  is  the  70  °C  data,  which  have  a  slope  larger  than 
the  data  at  60  °C.  Similarly,  in  Fig.  12  the  fits  to  Rreg  for  the 
average  of  all  the  cells  are  shown  for  the  four  test  tempera¬ 
tures.  Again,  the  fits  are  quite  good.  The  general  trend  is  the 
same  as  that  observed  for  Rdis  in  that  as  the  temperature  is 
increased,  the  slopes  and  intercepts  decrease.  The  exception 
is  the  70  °C  data.  Also,  note  from  Figs.  10-12  that  the  Rdis 
slopes  and  intercept  values  are  larger  than  those  for  Rreg  at 


the  same  test  temperature.  This  means  that  the  discharge 
resistance  is  increasing  at  a  greater  rate  than  the  regen 
resistance  at  the  same  temperature.  The  discharge  resistance 
at  zero  time  is  also  larger  than  the  zero  time  regen  resistance. 
This  type  of  analysis  was  done  on  each  of  the  cycle-life  test 
cells  at  the  A3,  A6,  A9%  SOC  swing  tests. 

A  functional  form  for  the  description  of  the  test  time, 
temperature,  SOC  and  A%  SOC  dependence  of  the  dis¬ 
charge  and  regen  resistances  was  assumed  to  be  essentially 
the  same  as  given  in  Eq.  (8)  for  the  calendar-life  tests: 

R(t,  T,  SOC,  A%  SOC))  =  A(T,  SOC,  A%  SOC)A2 

+  B(T,  SOC,  A%  SOC)  (8) 

with  the  addition  of  the  A  and  B  functions  being  dependent 
on  the  A%  SOC  swing  of  the  test  (either  A3,  A6  or  A9% 
SOC).  A(T. ;  SOC,  A%  SOC)  and  B(T,  SOC,  A%  SOC)  are 
assumed  to  be  only  functions  of  the  temperature,  the  SOC 
and  the  A%  SOC.  For  brevity,  A(T,  SOC,  A%  SOC)  is  simply 
referred  too  as  A,  and  B(J \  SOC,  A%  SOC)  is  referred  too  as 
B  unless  the  SOC  dependence  is  specifically  discussed.  Once 
the  functions  for  A  and  B  have  been  determined,  a  physical/ 
chemical  basis  for  the  functional  forms  are  attempted,  using 
the  fits  to  the  discharge  and  regen  resistances  as  a  guide. 

In  Fig.  13,  the  Rdis  and  Rreg  A  parameters  are  plotted  as  the 
exponential  of  the  inverse  of  the  temperature  (K),  i.e. 
A  =  a[exp(b/T )]  for  each  of  the  test  cells  (three  at  each  test 
temperature).  The  variation  in  the  actual  test  temperature, 
measured  by  a  thermocouple  mounted  on  the  cell  body,  is 
shown.  The  best-fit  fitting  parameters  are  shown  in  the  figure 
for  both  the  discharge  and  regen  A (7)  values.  From  the  fits 
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Fig.  11.  Average  of  all  INEEL  cycle-life  A3%  (60340,  60350,  60360  and  60370)  test  discharge  resistance  data  for  ATD  Gen  1  cells  plotted  as  a  function  of 
the  square-root  of  test  time.  The  model  fits  to  the  relation:  R(t,  T)  =  A(T)t1^2  +  B(T)  and  the  correlation  coefficient  ( R 2)  are  given  in  the  figure  for  each 
temperature. 
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Fig.  12.  Average  of  all  INEEL  cycle-life  A3%  (60340,  60350,  60360  and  60370)  test  regen  resistance  data  for  ATD  Gen  1  cells  plotted  as  a  function  of  the 
square-root  of  test  time.  The  model  fits  to  the  relation:  R(t,  T)  =  A(T)t1^2  +  B(T)  and  the  correlation  coefficient  ( R 2)  are  given  in  the  figure  for  each 
temperature. 


given  in  Fig.  13,  the  Rdis  A  parameter  is  seen  to  decrease  at  a 
slightly  greater  rate  with  increasing  temperature  than  does 
the  Rreg  A  parameter. 

Fig.  14  presents  a  similar  treatment  of  the  temperature 
dependence  of  the  B  parameter  as  determined  from  each  cell, 


also  for  the  (60340,  60350,  60360  and  60370)  test  condi¬ 
tions.  In  this  figure  are  plotted  the  Rdis  and  Rreg  B  parameters 
as  a  function  of  the  exponential  of  the  inverse  of  the 
temperature  (K),  i.e.  B  =  c[exp(d/7)].  The  actual  tempera¬ 
ture  of  the  cell  is  used  in  the  plots.  The  fitting  parameters  and 
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Fig.  13.  A(T)  parameters  from  data  fits  to  the  cycle-life  A3%  (60340,  60350,  60360  and  60370)  test  to  the  model  R(t,  T)  =  A(T)t F2  +  B(T)  are  shown  as  a 
function  of  the  inverse  test  temperature  (K).  Fitting  parameters  to  the  model  A(T)  =  a[exp(b/T )]  are  given  for  the  discharge  and  regen  resistance  data. 
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Fig.  14.  B(T)  parameters  from  data  fits  to  the  cycle-life  A3%  (60340,  60350,  60360  and  60370)  test  to  the  model  R(t ,  T)  —  A(7)A2  +  B(T)  are  shown  as  a 
function  of  the  inverse  test  temperature  (K).  Fitting  parameters  to  the  model  B(T)  =  c{exp(d/T)]  are  given  for  the  discharge  and  regen  resistance  data. 


R ~  are  given  in  the  figure.  Hence,  as  was  the  case  for  the  A 
parameter,  the  fits  to  the  temperature  dependence  of  the  B 
parameter  can  be  reasonably  correlated  by  the  exponential  of 
the  inverse  temperature  relation.  Note  that  the  scatter  in  the 
B  parameter  is  much  less  than  the  scatter  in  the  A  parameter. 
As  observed  in  this  study,  the  A  parameter  determined  from 
the  square-root  of  test  time  data  fits  is  very  sensitive  to  the 
data.  This  occurs  since  subtle  changes  in  the  slopes  of  the 
data  curves  fitted  to  the  square-root  of  the  test  time  can 
significantly  alter  the  A  parameter.  The  B  parameter,  which 
represents  the  intercept  of  the  square-root  of  test  time  fits,  is 
less  sensitive  to  the  test  data.  From  the  fits  given  in  Figs.  13 
and  14,  it  can  be  seen  that  the  Rdis  A  and  B  parameters 
decrease  at  a  slightly  greater  rate  with  increasing  tempera¬ 
ture  than  do  the  Rreg  A  and  B  parameters. 

Upon  close  examination  of  the  values  of  A  and  B  for 
the  discharge  and  regen  resistances  shown  in  Figs.  13  and  14, 
it  is  not  obvious  that  any  one  set  of  the  parameters  at  a 
given  temperature  should  be  eliminated.  Even  though  the 


discharge  and  regen  resistance  at  70  °C  appear  to  indicate 
that  a  new  mechanism  is  leading  to  a  resistance  increase  of 
the  cell  during  the  test  is  occurring,  all  of  the  resistance 
values  at  all  of  the  test  temperatures  were  used  in  the  model. 
The  values  for  the  discharge  and  regen  activation  energies 
and  the  pre-exponential  coefficients  determined  from  the 
data  fits  shown  in  Figs.  13  and  14  for  the  A3%  SOC  swing 
test  are  given  in  Tables  9  and  10.  Note  that  the  parameters  a , 
b ,  c  and  d  can  be  SOC-dependent  as  well  as  A%  SOC- 
dependent,  as  explicitly  shown  in  the  following  equation: 


R(t ,  T,  SOC,  A%  SOC)  =  a( SOC,  A%  SOC) 

~b(  SOC,  A%  SOC)' 


x  <  exp 


T 


4/2 


+  c(SOC,  A%  SOC) 

d(SOC,  A%  SOC) 


x  <  exp 


T 


(9) 


Table  9 

Values  of  Eact  at  60%  SOC  from  analysis  of  cycle-life  A3,  A6  and  A9%  SOC  test  data  using  the  relationship:  R(t,  T)  =  A(T)tl/2  +  B(T),  where  A(T)  = 
a[exp(EactA/RT)]  and  B(T)  =  c[exp(Eact?B/RT)] 


Test  condition 

Activation  energy  (Eac{A,  kJ/mol) 

Activation  energy  (Eact  B,  kJ/mol) 

Rdis  (60340,  60350,  60360  and  60370) 

18.2  (4.35) 

5.96  (1.42) 

tfreg  (60340,  60350,  60360  and  60370) 

13.8  (3.30) 

5.76  (1.38) 

Rdis  (60640,  60650,  60660  and  60670) 

30.2  (7.23) 

4.52  (1.08) 

Rreg  (60640,  60650,  60660  and  60670) 

19.0  (4.53) 

4.88  (1.17) 

RdiS  (60940,  60950,  60960  and  60970) 

11.0  (2.63) 

5.57  (1.33) 

tfreg  (60940,  60950,  60960  and  60970) 

9.51  (2.27) 

4.88  (1.17) 

Eact  are  activation  energies,  and  a  and  c  are  pre-exponential  constants.  The  values  given  in  parenthesis  are  in  kcal/mol. 
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Table  10 


Values  of  the  pre-exponential  constants  a  and  c  at  60%  SOC  from  analysis  of  cycle-life  test  data  analysis  at  A3,  A6%  and  A9%  SOC  using  the  relationship: 
R(t,  T)  =  A(T)tl/ 2  +  B(T),  where  A(T)  =  a[exp(EactiA/7?7)]  and  B(T)  =  c[Qxp(EactB/RT)] 


Test  condition 

1  /r\ 

Pre-exponential  constant  a  (Q  per  (day)  ) 

o 

Pre-exponential  constant  c  (x  10  Q) 

Rdis  (60340,  60350,  60360  and  60370) 

7.53  x  1CT7 

4.97 

Rreg  (60340,  60350,  60360  and  60370) 

2.14  x  1(T6 

3.72 

Rdis  (60640,  60650,  60660  and  60670) 

1.14  x  1(T8 

8.30 

Rreg  (60640,  60650,  60660  and  60670) 

3.65  x  1CT7 

5.22 

Rdis  (60940,  60950,  60960  and  60970) 

1.21  x  1CT5 

5.55 

Rreg  (60940,  60950,  60960  and  60970) 

1.25  x  10“5 

5.08 

Since  only  one  SOC,  60%,  was  studied  in  this  paper,  the 
SOC  dependence  of  the  various  parameters  could  not  be 
determined.  The  dependence  of  the  parameters  on  A%  SOC 
is  discussed  later  in  this  paper. 

The  comparison  between  the  model  predictions  using 
Eq.  (9)  with  the  parameter  values  as  given  in  Tables  9 
and  10  and  the  experimental  discharge  and  regen  resistance 
data  at  the  60%  SOC,  A3%  SOC  charge  swing  cycle-life 
test  are  presented  in  Figs.  15  and  16.  The  model  predictions 
are  quite  good  for  the  40  and  50  °C  data.  The  model, 
however,  would  predict  that  the  60  °C  resistance  would 
be  larger  than  the  70  °C  resistance.  Experimentally,  this  is 
not  the  case,  as  the  70  °C  resistance  values  are  slightly 
higher  than  the  60  °C  values.  Based  on  the  model  and  the 
experimental  data,  new  physical/chemical  processes  may  be 
occurring  at  60  and/or  70  °C  that  cause  the  resistance  to  not 
follow  the  general  trend  of  a  lower  resistance  the  higher  the 
temperature.  Also  shown  in  Figs.  15  and  16  are  the  model’s 


prediction  of  the  discharge  and  regen  resistances  for  a  test 
temperature  of  25  °C.  The  ability  to  predict  the  discharge 
and  regen  resistances  at  temperatures  other  than  ones  studied 
experimentally  is  one  of  the  main  points  in  conducting  this 
modeling  effort.  Unfortunately,  there  were  no  data  acquired 
at  25  °C  or  other  temperatures,  which  may  be  used  to 
validate  the  model,  particularly  at  the  lower  temperatures. 

The  discharge  and  regen  resistances  as  a  function  of  test 
time  at  temperature  for  the  A6  and  A9%  SOC  swing  tests 
were  analyzed  in  the  same  manner  as  was  used  for  the  A3% 
SOC  swing  test  described  earlier  [6].  The  best  fit  to  the  time 
dependence  of  the  increase  in  the  discharge  and  regen 
resistance  was  determined  to  be  a  square-root  of  time 
dependence  as  is  shown  in  Eqs.  (8)  and  (9).  The  A  and  B 
parameters  were  also  found  to  be  best  represented  by  an 
Arrhenius-type  of  temperature  function  as  shown  in  Eq.  (9). 
The  fits  of  the  model  as  given  in  Eq.  (9)  used  the  resistance 
values  at  all  of  the  test  temperatures  since  it  was  not  obvious 
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Fig.  15.  Model  predictions  compared  to  cycle-life  A3%  discharge  resistance  data  for  ATD  Gen  1  (60340,  60350,  60360  and  60370)  tests.  Model  used  the  a,  b, 
c  and  d  parameters  determined  when  all  the  test  temperatures  were  used  to  fit  the  resistance  data.  The  model  prediction  for  a  test  temperature  of  25  °C  is  also 
shown. 
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Fig.  16.  Model  predictions  compared  to  cycle-life  A3%  regen  resistance  data  for  ATD  Gen  1  (60340,  60350,  60360  and  60370)  tests.  Model  used  the  a,  b,  c 
and  d  parameters  determined  when  all  the  test  temperatures  were  used  to  fit  the  resistance  data.  The  model  prediction  for  a  test  temperature  of  25  °C  is  also 
shown. 


which  temperature  data  to  eliminate  from  the  plots  of  the 
A (T)  and  B(T)  functions  as  a  function  of  the  inverse  of  the 
test  temperature.  The  values  of  the  discharge  and  regen 
activation  energies  and  the  pre-exponential  coefficients 
determined  from  this  analysis  are  given  in  Tables  9  and 
10.  The  average  of  the  experimental  discharge  and  regen 
resistance  and  the  model  fits  to  the  data  using  Eq.  (9)  with 
the  parameter  values  given  in  Tables  9  and  10  are  shown  in 
Figs.  17  and  18  for  the  A6%  SOC  swing  test,  and  in  Figs.  19 
and  20  for  the  A9%  SOC  swing  test,  respectively.  Also 
shown  are  the  model  predictions  for  the  discharge  and  regen 
resistance  if  the  tests  were  conducted  at  25  °C.  The  model 
fits  are  fairly  good  for  the  40  and  50  °C  discharge  and  regen 
resistance.  The  correlation  to  the  60  and  70  °C  data  is  not  as 
good,  but  the  model  predictions  do  bracket  the  resistance  at 
these  temperature  values. 

Upon  examination  of  the  activation  energies  presented  in 
Table  9,  several  overall  trends  can  be  observed.  The  dis¬ 
charge  activation  energies  for  Ea cuA  are  greater  than  the 
corresponding  regen  EactA  activation  energies.  From  Eq.  (9), 
this  indicates  that  the  time  rate  of  change  of  the  discharge 
resistance  is  greater  than  the  regen  resistance.  There  appears 
to  be  no  general  trend  for  the  EactA  activation  energies, 
except  to  observe  that  the  A6%  SOC  test  values  are  higher 
than  those  for  A3  and  A9%  test  values.  The  Eact  B  activation 
energies  are  less  than  the  EactA  activation  energies.  The 
EactB  activation  energies  are  all  about  the  same,  probably 
within  experimental  error.  At  present,  the  authors  know  of  no 
reason  why  the  values  given  in  Table  9  are  what  they  are. 
They  are  about  the  correct  order  of  magnitude  for  processes 


that  may  occur  in  a  battery  such  as  the  various  diffusion 
processes  that  appear  to  have  activation  energies  of  several 
tens  of  kJ/mol.  The  authors  have  been  unable  to  find 
literature  values  for  lithium-ion  diffusion  through  the  SEI 
layer  present  on  the  anode  or  cathode,  or  its  diffusion  in  the 
anode  or  cathode  materials  for  the  electrode  and  electrolyte 
compositions  used  in  the  ATD  Gen  1  cells.  Similarly,  the 
pre-exponential  constant  a  shown  in  Table  10  is  smaller  for 
the  discharge  resistance  than  the  regen  resistance.  The  pre¬ 
exponential  constant  c  for  the  discharge  resistance  is  some¬ 
what  greater  than  the  regen  resistance  value  at  each  A%  SOC 
test.  Why  this  is  so  is  not  presently  known. 

In  Eq.  (9),  it  was  assumed  that  the  resistance  was  a 
function  of  the  SOC  of  the  cell  at  the  beginning  of  the  test 
pulse  sequence,  and  that  it  was  also  a  function  of  the  SOC 
swing  during  the  test,  A%  SOC.  As  only  one  SOC  was 
tested,  60%  SOC,  the  dependence  of  the  resistance  on  this 
variable  could  not  be  determined.  Attempts  were  made  to 
find  a  consistent  model  of  the  parameters  a ,  b ,  c  and  d  as  a 
function  of  the  A%  SOC  of  the  test.  No  reasonably  consistent 
model  was  found,  nor  were  any  insights  into  the  A%  SOC 
dependence  found  in  the  literature.  Physical/chemical  diag¬ 
nostics  being  conducted  on  the  test  cells,  as  part  of  the  ATD 
program,  may  provide  future  insight  into  the  processes 
responsible  for  these  model-determined  values. 

3.4.  Comparison  between  calendar-  and  cycle-life  tests 

Upon  comparison  of  the  measured  discharge  and  regen 
resistance  values  for  the  calendar-life  tests  conducted  at  60 
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Fig.  17.  Model  predictions  compared  to  cycle-life  A6%  discharge  resistance  data  for  ATD  Gen  1  (60640,  60650,  60660  and  60670)  tests.  Model  used  the  a,  b, 
c  and  d  parameters  determined  when  all  the  test  temperatures  were  used  to  fit  the  resistance  data.  The  model  prediction  for  a  test  temperature  of  25  °C  is  also 
shown. 


and  80%  SOC  with  those  measured  at  60%  SOC  for  the 
cycle-life  tests  (Figs.  5,  7-9  and  15-20),  it  can  be  seen  that 
the  both  the  discharge  and  regen  resistance  were  greatest  for 
the  80%  SOC  calendar-life  test.  The  discharge  and  regen 
60%  SOC  calendar-life  test  resistances  were  also  greater 


than  the  60%  SOC  cycle-life  tests  at  A3,  A6  and  A9%  at 
times  less  than  4  weeks.  This  observation  is  counter-intuitive 
as  one  would  had  expected  that  the  cycle-life  tests  at  A3,  A6 
and  A9%  SOC  would  have  been  more  “abusive”  on  the  cells 
resulting  in  a  greater  increase  in  their  resistance  with  test 
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Fig.  18.  Model  predictions  compared  to  cycle-life  A6%  regen  resistance  data  for  ATD  Gen  1  (60640,  60650,  60660  and  60670)  tests.  Model  used  the  a,  b,  c 
and  d  parameters  determined  when  all  the  test  temperatures  were  used  to  fit  the  resistance  data.  The  model  prediction  for  a  test  temperature  of  25  °C  is  also 
shown. 
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Fig.  19.  Model  predictions  compared  to  cycle-life  A9%  discharge  resistance  data  for  ATD  Gen  1  (60940,  60950,  60960  and  60970)  tests.  Model  used  the  a ,  b, 
c  and  d  parameters  determined  when  all  the  test  temperatures  were  used  to  fit  the  resistance  data.  The  model  prediction  for  a  test  temperature  of  25  °C  is  also 
shown. 


time.  The  calendar-life  test  had  a  percentage  SOC  swing  of 
1%  (Table  1)  on  discharge  and  was  conducted  only  once  per 
day  while  the  cycle-life  tests  were  conducted  almost  con¬ 
tinuously  and  had  percentage  SOC  swings  of  3,  6  and  9%.  It 


was  also  initially  thought  that  the  order  of  the  resistance 
increase  for  the  cycle-life  tests  would  have  been  A9% 
SOC  >  A6%  SOC  >  A3%  SOC  because  of  the  greater 
perturbation  of  the  cell  the  greater  the  A%  SOC  swing.  This 
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Fig.  20.  Model  predictions  compared  to  cycle-life  A9%  regen  resistance  data  for  ATD  Gen  1  (60940,  60950,  60960  and  60970)  tests.  Model  used  the  a,  b,  c 
and  d  parameters  determined  when  all  the  test  temperatures  were  used  to  fit  the  resistance  data.  The  model  prediction  for  a  test  temperature  of  25  °C  is  also 
shown. 
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also  was  not  the  case  as  the  resistance  at  A3%  was  observed 
to  increase  to  a  greater  extent  that  at  A6  and  A9%  (at  least  for 
test  times  up  to  200  h  as  shown  in  Figs.  15-20).  It  may  well 
be  that  if  the  formation  of  the  SEI  layer  is  the  primary  cause 
of  the  resistance  increase,  that  the  rate  of  its  formation  is 
fairly  slow  such  that  the  cycle-life  tests  are  disrupting  its 
formation  on  a  rather  continuous  basis  due  to  the  dischar¬ 
ging/charging  cycles  compared  to  when  the  cell  is  sitting 
for  a  longer  period  of  time  between  discharging/charging 
cycles  as  in  the  calendar-life  tests.  That  this  may  be  the  case 
is  supported  by  the  pre-exponential  constant  a  for  the  tests 
as  are  given  in  Tables  8  and  10  for  the  various  tests.  It  is 
observed  in  these  tables  that  the  pre-exponential  constant  a 
for  the  60%  SOC  calendar-life  test  are  greater  than  for  the 
60%  SOC  cycle-life  tests.  The  time-rate-of-change  of  the 
discharge  and  regen  resistance  for  the  models  given  in 
Eqs.  (7)  and  (9)  is: 


d[R(t,  T,  SOC,  A%  SOC)] 
d  t 


a(  SOC,  A%  SOC) 

f  exp[/?(SOC,  A%  SOC) /T\  \ 

X  \  2tW  / 

(10) 


This  equation  shows  that  the  time-rate-of-change  of  the 
resistance  not  only  depends  on  the  a( SOC,  A%  SOC)  pre¬ 
exponential  constant,  but  also  on  the  b( SOC,  A%  SOC) 
parameter  that  is  related  to  the  activation  energy  of  the 
process.  The  activation  energies,  as  shown  in  Tables  7  and  9, 
are  greater  for  the  cycle-life  tests  than  for  the  calendar- 
life  tests,  so  the  actual  rate  of  increase  of  the  discharge  and 


regen  resistance  depend  on  both  the  a{ SOC,  A%  SOC)  pre¬ 
exponential  constant  and  also  on  the  value  of  the  term 
exp[£act(SOC,  A%  SO C)/RT].  This  relation  indicates  that 
an  increase  in  one  parameter  can  be  compensated  for  by  a 
decrease  in  the  other.  Direct  comparisons  of  the  discharge 
and  regen  resistance  for  the  80  and  60%  SOC  calendar-life 
tests  and  the  60%  SOC  cycle-life  tests  as  predicted  using  the 
model  as  given  in  Eqs.  (7)  and  (9),  and  the  fitting  parameters 
given  in  Tables  5-10  for  a  test  temperature  of  25  °C  are 
given  in  Figs.  20  and  21,  respectively,  for  a  time  period  of 
365  days.  Examining  the  discharge  resistance  for  the 
(80C25)  test,  it  can  be  seen  that  the  it  is  the  greatest  at  this 
test  condition  compared  to  any  of  the  other  tests.  The 
discharge  resistance  for  the  (60C25)  test  condition  is  also 
greater  than  the  all  of  the  A%  SOC  cycle-life  tests  except  for 
the  A6%  test  whose  resistance  is  predicted  to  becomes  larger 
after  ~100  days.  The  A3%  SOC  cycle-life  discharge  resis¬ 
tance  is  greater  than  the  A9%  SOC  tests  over  the  entire  time 
period  shown  in  Fig.  21.  The  model  predictions  for  the  regen 
resistances  at  25  °C  (Fig.  22)  for  the  various  tests  shows  that 
the  calendar-life  test  caused  the  resistance  to  increase  to  a 
greater  extent  then  any  of  the  A%  SOC  tests.  At  the  longer 
test  times  (~>15  days),  the  regen  resistance  for  the  cycle-life 
tests  is  in  the  order  A6%  SOC  >  A3%  SOC  >  A9%  SOC.  At 
present,  we  have  no  explanation  for  these  trends  as  they  are 
not  what  would  have  been  expected  before  the  actual  testing 
of  the  batteries.  Thus,  the  exact  nature  of  the  test,  and  the 
SOC  at  which  the  test  is  conducted  (and  the  A%  SOC  swing 
for  the  cycle-life  tests)  all  have  a  major  influence  on  the 
aging  of  the  batteries  and  the  change  in  the  discharge  and 
regen  resistance  during  a  given  test. 
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Fig.  21.  Model  predictions  for  a  time  period  of  365  days  for  calendar-  and  cycle-life  discharge  resistance  for  ATD  Gen  1  cells  at  different  test  conditions. 
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Fig.  22.  Model  predictions  for  a  time  period  of  365  days  for  calendar-  and  cycle-life  regen  resistance  for  ATD  Gen  1  cells  at  different  test  conditions. 


4.  Conclusions 

This  paper  presented  the  test  results  pertaining  to  a  group 
of  ATD  program  Gen  1  lithium-ion  batteries.  The  cells 
underwent  a  number  of  electrical  performance  tests  to 
determine  their  electrochemical  performance  at  25  °C.  Spe¬ 
cial  calendar-  and  cycle-life  tests  were  also  conducted  at 
elevated  temperatures  of  40,  50,  60  and  70  °C  for  specified 
periods  of  time. 

One  of  the  specific  tests  for  which  the  data  are  presented 
and  discussed  was  a  special  calendar-life  test  conducted 
once  per  day  for  specified  periods  of  time,  depending  on  the 
test  temperature.  The  calendar-life  test  was  conducted  at  60 
and  80%  SOC.  During  the  calendar-life  test,  the  discharge 
resistance  was  determined  from  the  discharge  portion  of  the 
test  using  Ohm’s  law.  The  regen  resistance  was  determined 
from  the  regen  portion  of  the  test  also  using  Ohm’s  law. 

General  observations  derived  from  this  study  are  the 
following. 

1.  Both  the  discharge  and  regen  resistances  have  a  non¬ 
linear  increase  with  respect  to  time  at  test  temperature. 

2.  The  discharge  resistance  is  greater  than  the  regen 
resistance  at  all  test  temperatures  of  40,  50,  60  and 
70  °C. 

3.  For  both  the  discharge  and  regen  resistance,  generally 
the  higher  the  test  temperature,  the  lower  the  resistance. 

4.  The  70  °C  discharge  and  regen  resistance  data  at  80% 
SOC  do  not  follow  the  general  trend  of  the  rest  of  the 
data  in  that  the  resistance  at  this  temperature  is  slightly 
greater  than  at  60  °C.  At  60%  SOC,  the  discharge  and 


regen  resistance  data  indicate  a  process  is  occurring  that 
causes  the  60  °C  resistance  to  be  slightly  greater  or 
comparable  to  those  at  50  °C.  At  this  SOC,  the  70  °C 
may  also  have  been  influenced.  These  observations 
appear  to  indicate  that  a  new  physical/chemical  process 
is  occurring  that  causes  the  resistance  to  increase  in  an 
anomalous  manner.  They  also  indicate  that  the  SOC  at 
which  the  test  was  conducted  may  influence  the 
temperature  at  which  the  onset  of  these  new  processes 
occurs.  The  exact  nature  of  these  processes  is  not 
presently  known. 

5.  Both  the  discharge  and  regen  resistances  at  a  given 
temperature  are  greater  at  80%  SOC  than  they  are  at 
60%  SOC. 

6.  There  are  also  differences  in  the  rate  of  increase  of  the 
resistances  in  that  the  80%  SOC  resistance  increases 
with  time  at  test  temperature  faster  than  at  60%  SOC. 

A  life  model  was  developed  to  account  for  the  time,  tem¬ 
perature  and  SOC  of  the  batteries  during  the  calendar-life  test. 
The  of  modeling  effort  led  to  the  relationship  for  the  discharge 
and  regen  resistance  having  the  form:  R(t ,  7\  SOC)  = 
atSOCKexpt^SOCyr]}/1/2 3 4  +  c(SOC){exp[J(SOC)/7]}.  The 
square-root  of  time  dependence  can  be  accounted  for  by  either  a 
one-dimensional  diffusion-type  of  mechanism,  presumably  of 
the  lithium-ions  or  by  a  parabolic  growth  mechanism  for  the 
growth  of  a  thin-film  SEI  layer  on  the  anode  and/or  cathode. 
The  temperature  dependence  of  the  resistances  have  the  form  of 
an  Arrhenius-like  process  that  could  be  related  to  the  diffusion 
of  the  lithium-ions  or  to  chemical  processes  related  to  the 
formation  of  the  SEI  surface  film.  The  function  R(t,  T,  SOC) 
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was  then  be  used  to  correlate  the  discharge  and  regen  resistance 
data  and  to  predict  what  the  resistances  would  be  at  different 
test  times  at  particular  test  conditions  and  at  different  test 
temperatures. 

The  other  specific  test  for  which  the  data  are  presented  and 
discussed  in  this  paper  was  a  special  cycle-life  test  con¬ 
ducted  for  a  4- week  period  at  test  temperatures  of  40,  50  and 
60  °C,  and  for  a  2- week  period  at  70  °C.  This  test,  consisting 
of  three  specified  discharge  and  charge  protocols,  was 
designed  to  have  a  A3,  A6  or  A9%  cumulative  SOC  swings 
during  the  discharge  portion  of  the  test  cycle. 

The  results  of  the  cycle-life  testing  indicate  that  both  the 
discharge  and  regen  resistance  increased  in  a  non-linear 
manner  as  a  function  of  the  test  time  at  each  A%  SOC  test. 
The  magnitude  of  the  discharge  and  regen  resistances 
depended  on  the  temperature  and  A%  SOC  at  which  the 
test  was  conducted.  General  observations  derived  from  this 
study  are  the  following. 

1.  Both  the  discharge  and  regen  resistances  have  a  non¬ 
linear  increase  with  respect  to  time  at  test  temperature, 
i.e.  as  the  number  of  test  cycles  increased,  the  discharge 
and  regen  resistance  increased  also. 

2.  For  a  given  A%  SOC  test,  the  discharge  resistances  are 
greater  than  the  regen  resistance  at  all  of  the  test 
temperatures  of  40,  50,  60  and  70  °C. 

3.  For  both  the  discharge  and  regen  resistances,  generally 
the  higher  the  test  temperature,  the  lower  the  resistance. 

4.  The  70  °C  discharge  and  regen  resistance  data  do  not 
always  follow  the  general  trend  of  the  rest  of  the  data  in 
that  the  resistance  at  this  temperature  are  slightly  greater 
than  at  60  °C.  This  observation  appears  to  indicate  that 
new  physical/chemical  processes  are  occurring  that 
causes  the  resistance  to  increase.  The  exact  nature  of 
these  processes  is  not  presently  known. 

5.  At  each  of  the  four  test  temperatures,  the  magnitude  of 
the  discharge  and  regen  resistance  was  in  the  following 
order  for  the  first  200  h  test  period:  A3%  SOC  >  A9% 
SOC  >  A6%  SOC.  At  longer  test  times,  the  model 
predicts  that  the  relative  order  of  the  discharge  and  regen 
resistance  at  a  given  A%  SOC  test  would  be  A6% 
SOC  >  A3%  SOC  >  A9%  SOC.  No  explanation  is 
currently  known  for  this  observation.  There  are  research 
groups  within  the  ATD  program  currently  conducting 
physical/chemical  studies  on  the  test  cells  that  may  offer 
insight  into  this  observation. 

A  life  model  was  developed  to  account  for  the  time, 
temperature,  SOC  and  A%  SOC  of  the  batteries  during 
the  cycle-life  test.  The  functional  form  of  the  model  is  given 
by:  R(t ,  T,  SOC,  A%  SOC)  =a(SOC,  A%  SOC){exp[Z?(SOC, 
A%  SOC)/:T]}r1/2  +  c(SOC,  A%  SOC){exp[d(SOC,  A% 
SOC)/7]}.  This  model  is  essentially  the  same  as  was  used 
to  model  the  calendar-life  discharge  and  regen  resistances. 
The  values  determined  for  the  functions  a,  b,  c  and  d  were 
different  than  those  determined  for  the  calendar-life  tests. 
Due  to  the  lack  of  testing  at  SOC  values  other  than  60% 


SOC,  the  exact  form  of  the  SOC  dependence  could  not  be 
determined  from  the  cycle-life  data.  Attempts  were  made  to 
find  a  consistent  correlation  of  the  observed  resistance 
changes  with  the  A%  SOC  of  the  tests,  but  no  model  based 
on  physical/chemical  processes  was  found. 
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